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Abstract. We consider the finite-temperature frequency and momentum dependent 
two-point functions of local operators in integrable quantum field theories. We focus 
on the case where the zero temperature correlation function is dominated by a delta- 
function line arising from the coherent propagation of single particle modes. Our 
specific examples are the two-point function of spin fields in the disordered phase of 
the quantum Ising and the 0(3) nonlinear sigma models. We employ a Lehmann 
representation in terms of the known exact zero-temperature form factors to carry 
out a low-temperature expansion of two-point functions. We present two different but 
equivalent methods of regularizing the divergences present in the Lehmann expansion: 
one directly regulates the integral expressions of the squares of matrix elements in 
the infinite volume whereas the other operates through subtracting divergences in a 
large, finite volume. Our central results are that the temperature broadening of the line 
shape exhibits a pronounced asymmetry and a shift of the maximum upwards in energy 
("temperature dependent gap"). The field theory results presented here describe the 
scaling limits of the dynamical structure factor in the quantum Ising and integer spin 
Heisenberg chains. We discuss the relevance of our results for the analysis of inelastic 
neutron scattering experiments on gapped spin chain systems such as CsNiCl3 and 
YBaNi0 5 . 
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1. Introduction 

Progress over the last two decades [1, 2, 3, 4, 5, 6, 7, 8, 9, 10] has made it possible to 
determine zero temperature dynamical response functions in massive integrable models 
of quantum field theory (QFT) by means of the "form factor bootstrap approach". 
More precisely, response functions can be calculated exactly at low frequencies (several 
times the mass gap) and to high accuracy at intermediate frequencies. The T = 
dynamics described by integrable QFTs exhibits a number of interesting phenomena 
such as dynamical mass generation, spin-charge separation, and other kinds of quantum 
number fractionalization. The results obtained by the form factor bootstrap approach 
have important applications in condensed matter systems [11] such as quantum magnets 
[12], Mott insulators [13], doped ladder materials [14], carbon nanotubes [15] and ultra- 
cold atomic gases [16]. 

At the heart of the form factor bootstrap approach is the notion that in an 
integrable model the scattering of elementary excitations is purely elastic by virtue 
of the existence of an infinite number of local conservation laws. There is no particle 
production and the individual particle momenta are conserved in scattering events. 
Correspondingly, the scattering of n particles in an integrable model can always be 
reduced to a sum of two-body scattering events. The resulting simplified nature of 
the exact Hamiltonian eigenstates in a massive integrable field theory permits the 
computation of zero temperature correlation functions as follows. One first employs 
a Lehmann representation in terms of n-particle Hamiltonian eigenstates \n, {s n }) with 
energy i£[{s n }], where {s n } labels the corresponding sets of good quantum numbers 

X o(r,x) =(0|T T O(r,x)O t (0,0)|0) 
r=o 

oo 

= 5^e- Tjs M->]<0|O(x,0)|n;K}>(n;K}|Ot(0,0)|0> (r > 0). (1) 

n=0 

As the energies are simply given as sums over the single particle energies of the 
elementary excitations, we have reduced the computation of the response function to 
computing a set of matrix elements or "form factors", (0\O(t, x)\n; {s n }). 

To perform this computation we again exploit integrability. The ability to express 
\n, {s n }) as a collection of n distinct particles allows one to write down a set of algebraic 
constraints that the matrix elements must satisfy [1] . These constraints encode both the 
simplified form the scattering of n particles takes in an integrable model together with 
analytic constraints coming from crossing symmetries present in a relativistic quantum 
field theory. While these constraints can be written down for eigenstates involving 
an arbitrary number, n, of particles, they become increasingly cumbersome to solve 
as n increases. Fortunately, if we are interested merely in the behavior of the low 
energy spectral function, it is only necessary to compute matrix elements involving 
a few particles. The spectral function is obtained by Fourier transforming (1) with 
respect to space and imaginary time and then analytically continuing to real frequencies 
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iu n — > uj + iO) 

-. oo 

-Im X o(u, q) = ^J2Jl{\ M) I' 6 ^ ~ E ti s ^V S(q - P[{s n }]) 



7T 

n=0 



- 6 I (n; {s n }\O^(0, 0)|0) I 2 % + E Sn ) 6(q + P[{a B }]) }, (2) 



where -P[{s n }] denotes the momentum of |n, {s n }) and e = ±1 depending on whether 
the field O is bosonic or fermionic. The presence of the delta function in the above 
expression for the spectral function guarantees that at an energy u only eigenf unctions 
with this exact energy contribute. As the theory is massive, eigenstates with n-particles 
will have a minimum energy, nA (supposing the particles all have mass A). Thus for 
energies, u> < nA, such eigenstates will not contribute to the spectral function. For low 
energies, only a small finite number of matrix elements need to be computed in order 
to obtain exact results for the spectral function. Even at higher energies, it has been 
typically found that the sum of matrix elements is strongly convergent, and that matrix 
elements with higher particle number make only an extremely small contribution to the 
spectral function [3, 17, 18, 19, 11]. 

While the above described approach has been successful at computing zero 
temperature correlation functions, it is not a settled question whether the form factor 
bootstrap approach can be used universally to gain information about T > dynamical 
correlations. This question has been investigated in a number of specific instances 
[20, 21, 22, 23, 24]. On the basis of this past work, there appear, in various forms, two 
main difficulties in doing so. These difficulties become particularly acute in theories 
which are interacting. To illustrate these problems, we write out the corresponding 
form-factor expansion for a Green's function at finite temperature: 

--Im Xo (w, 9) = ^^- E ^ S ^ + ^K*™}]) % ~ P H S n}} + P[{Sm}]) 

ns n 
msra 

x [ e -PE[{sm}) _ ee -mW}]] | (m; { Sm }\O(0,0)\n; K}>| 2 . (3) 

The first difficulty can be seen in that the form-factor expansion now involves two 
sums over eigenstates, the first sum arising from the insertion of the resolution of the 
identity as before and the second sum coming from the Boltzmann trace associated 
with working at finite temperature. Concomitantly, working at a particular energy, 
u, no longer guarantees that only a finite number of matrix elements will contribute 
to the spectral function. This problem was partially resolved in Ref. [23]. There it 
was advocated that the Boltzmann factor provides a natural small parameter, that is, 
terms involving eigenstates \m; {s m }) with many particles (and so i?[{s m }] large) play 
only a specific, limited role, at least at low temperatures. It was shown there that for 
a certain class of correlation functions, it was possible to develop a low temperature 
expansion of the response function. However this work left unresolved how to perform 
the low temperature expansion in general. In particular, it did not address how to 
develop the low temperature expansion for response functions with singular features at 
zero temperature. 
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Such response functions form a class of great physical interest. As an important 
physical example, the spin response, S(u,q), of a gapped quantum spin chain as 
represented by either the quantum Ising model or the 0(3) NLSM is given by S(u, q) = 
Z(q)8(uj — e(q)) + . . . at zero temperature, where the dots indicate very weak multi- 
paticle scattering continua. At finite temperature this (^-function response broadens out 
to a lineshape with finite width. But it has been unclear how to capture this broadening 
in the context of the form factor expansion as given in (3). It is one of the achievements 
of this paper to demonstrate how this can be done. 

The second difficulty is that unlike the zero temperature case, we are now faced 
with computing matrix elements between states, |n; {s n }) and |m; {s m }) both with finite 
particle number (as opposed to amplitudes (0|C(0, 0)|n; {s n }) governing the transition 
between the vacuum and some eigenstate). In an infinite volume, matrix elements 
of the form (m; {s m }\O(0, 0)|n; {s n }) with n, m ^ possess singularities. Moreover 
these singularities are 'squared' as it is the quantity \(m; {s m }\0 (0,0) \n; {s n })\ 2 that 
determines the correlation function. Making sense of these squared singularities amounts 
to giving a sensible interpretation to terms of the form, "5(0)". 

It is perhaps the main achievement of this paper that we have demonstrated a set of 
regularization procedures for such singularities. These ultimately arise as a consequence 
of working in an infinite volume, where the momenta of particles found in the n-particle 
and m-particle states, |n; {s n }) and \m; {s m }), can be identical. In a matrix element 
of the form, (m; {s m }\0 (0,0) \n; {s n }), identical momenta lead to divergences. One 
approach is thus to work in a large but finite volume. In such cases the momenta 
(at least in the relevant examples) of the particles composing the states \n; {s n }) and 
\ m ]{s m }) are never equal leaving the matrix element (m; {s m }\0 (0,0) \n; {s n }) finite. 
Such an approach is feasible as (m; {s m }|C(0, 0)|n; {s n }) retains the same momentum 
dependencies as in an infinite volume. Provided we work in a large, finite volume, 
the sole difference in evaluating the form-factor (m; {s m }\0 (0,0) \n; {s n }) for general 
momenta between the finite and infinite volume cases lies in taking into account that 
in finite volume the momenta of the states |n; {s n }) and \m; {s m }) are quantized [25]. 
Thus we can still use in a large, finite volume the infinite volume constraints that govern 
these matrix elements. 

The use of finite volume regularization has precedent. For the special case of the 
quantum Ising model this particular problem has been solved by calculating the form 
factors on the cylinder [26, 27] - that is, the form factors are computed in an arbitrary, 
not merely asymptotically large, finite volume. These results have then been used in 
Ref. [24] to recover the semiclassical expression obtained previously by Sachdev and 
Young [28] and in Ref. [29, 30] to carry out low and high temperature expansions 
for the spin-spin correlation function. Form factors as computed on the cylinder have 
also been exploited to derive partial differential equations governing finite temperature 
correlation functions in the quantum Ising model [31]. In this work, the large space and 
time asymptotics of spin correlations were determined. Such partial different equations 
for correlation functions have been derived in a number of instances [32] and employed 
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very recently [33] to analyze time dependent zero temperature correlators in the Ising 
model. However like in Ref. [31], the underlying theory always has had trivial scattering, 
that is the scattering matrix is momentum independent and diagonal. It is presently 
not known how to generalize the results of [26, 27] to integrable QFTs with non-trivial 
S-matrices. 

In addition to our use of a finite volume regularization scheme, we demonstrate in 
this paper a new regularization technique that operates in an infinite volume. In this 
regularization scheme, ambiguous terms such as "5(0)" are absent by construction. We 
demonstrate explicitly that our two regularization schemes lead to the same answers. 
We do so in two examples, the quantum Ising model and the 0(3) non-linear sigma 
model (NLSM). That the latter model is interacting and so highly non-trivial provides 
a strong indication that our infinite volume regularization scheme is robust and so 
provides a candidate for a regulator that works generally. In particular it would be 
useful to check the scheme for cases involving form factors between two 2-particle states 
or a 2-particle and a 3-particle state. 

Our work, both in how we develop a low temperature expansion, and in our 
particular choice of integrable models, the quantum Ising model and the 0(3) NLSM, 
is motivated in large degree by recent inelastic neutron scattering experiments on 
several quantum magnets. A key objective of this recent experimental work has been 
to investigate how the spin dynamics crosses over from the strongly correlated zero 
temperature quantum regime to the classical high temperature regime [34, 35, 36]. In a 
system such as the spin-1 Heisenberg chain that supports a coherent, gapped, magnetic 
single-particle excitation at T = 0, the question arises of how the dominant feature in 
the dynamical structure factor, a delta- function at zero temperature, broadens at finite 
temperatures. As field theories, the quantum Ising model and the 0(3) NLSM describe 
the scaling limits of the (non-integrable) spin-chain Hamiltonians used to model these 
various experiments. 

Applying the methodology detailed in this paper, we analyze this finite temperature 
lineshape. As our central finding in this regard, we demonstrate that the lineshape 
is always asymmetric in energy, a feature that becomes more pronounced as the 
temperature increases. For the Ising model we further demonstrate the emergence of 
a "temperature dependent gap". A subset of our results on the lineshape have been 
previously reported in Ref. [37]. 

At temperatures T far below the spin gap A our technique complements previous 
semiclassical approaches to the study of the lineshape of quantum spin chains[28, 38, 39, 
40, 41]. Both in our methodology and in semi-classical approaches, the lineshape has 
been shown to be essentially Lorentzian for T C A and for energies, u, in the vicinity 
of the gap, A. However in our approach we can both study temperatures where the 
semi-classics is inaccurate as well as the entire lineshape, not merely uj ~ A. 

The outline of this paper is as follows. In Section 2 we set out our framework 
for deriving low-temperature expansions of dynamical correlation functions in massive 
integrable quantum field theories. In particular, we summarize how the space of 
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Hamiltonian eigenstates in an integrable model is handled as well as how we develop the 
low temperature expansion appropriate for computing the finite temperature lineshape. 
In Section 3 we apply the method to the quantum Ising chain. We introduce how 
we regularize the squares of form-factors both by working in finite volume as well as 
in our new infinite volume scheme. We show that in the context of the Ising model, 
these schemes are equivalent. In the following section, Section 4, we present a detailed 
discussion of the results obtained for the quantum Ising model as well as a comparison to 
the semiclassical results of Sachdev and Young. In the next two sections we move on to 
the 0(3) NLSM showing that our methodology also works for interacting theories with 
non-trivial (even non-diagonal) scattering matrices. In Section 5 we consider the case 
of the retarded Green's function of the vector field in the 0(3) nonlinear sigma model, 
the quantity that corresponds to the spin response of a gapped Heisenberg spin chain 
near wavevector tt. We again show that the two regularization schemes yield the same 
result for this correlation function. In Section 6 we detail the results that so arise for the 
low-temperature dynamics for the 0(3) NLSM and compare them to the semiclassical 
results of Ref [39]. The final section, Section 7, presents a summary and discussion 
of our results. Computational details on our new method for regularizing form factor 
squares directly in the infinite volume are presented in several appendices. 

2. General Formalism 

A defining feature of integrable quantum field theories is a basis of scattering states of 
"elementary" excitations, which are eigenstates of the Hamiltonian. It is customary to 
construct these states from the so-called Faddeev-Zamolodchikov algebra 



Here 8\ }2 are rapidity variables, a, b are quantum numbers and S is the exact two-particle 
scattering matrix describing the purely elastic scattering of the elementary excitation. 
The S-matrix is a solution to the Yang-Baxter equation, which can be thought of as 
a consistency condition for factorizable three-particle scattering. Using the Faddeev- 
Zamolodchikov operators, a Fock space of states can be constructed as follows. The 
vacuum is defined by 



Multiparticle states are then obtained by acting with strings of creation operators Z* b {ff) 
on the vacuum 



Z a (d x )Z h {Q 2 ) 

Za{ei)zl{B 2 ) 



Si b '(9 1 -9 2 )Z bl (e 2 )Z af (9 1 ) , 

si b '(e 1 -e 2 )zl(e 2 )zt,(e 1 ), 

2ir6 ab 6{0 1 - 9 2 ) + - 9 2 )Zl{6 2 )Z a ,{B 1 ). 



(4) 



z a (e)\o) = o . 



(5) 




(6) 



Energy and momentum of the states (6) are by construction additive 



s 



E S (9 1 ,...,9 S ) = J2< d i) ' = Acoshfl 
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P S (9 1 ,...,9 S ) = ^-sinh% ■ ( ? ) 

3=1 

In terms of this basis the resolution of the identity is given by 

dOx . . . M 



i = |o)(oi + EE / 

n=l {a;} J ~ ( 



(2ir) n n\ 



e n ...e 1 ) an ... a ^- a -{e 1 ...e n \ . (8) 



In the basis of scattering states introduced above, the following formal spectral 
representation for the retarded finite temperature two-point function of the local 
operator O holds 

xoMH^EO^)- (9) 

r,s=0 

Here C r>s denotes the contribution with r particles in the thermal trace and s in the 
intermediate state 

p/3 poo 

C? s (u,q)= / dr dxe iuJ " T - iqx C r , s (r,x) 



iu) n —*uj+iO 

(- ~\ - ST f d0 ^--- d0 r f d6 'l ■ ■ - de 's „-(3E r „-T(E a -E r ) 

rA ' ] ~ s \iJ (2tt)T! J (2nYs\ 

x e -*Vr- p ')* \ a '- a r(6 1 ...6 r \O{0,0)\6' s ...6' 1 ) a ,.. < \ 2 . (10) 
The partition function can formally be expressed as 

z = (o\o) + Y,[ t e ~" <e) b{mb 

b J n 

+ E / d ^%e'^<^ ^ W 2 + • • • = f>- (11) 

61,62 ^ ' n=0 

Both the partition function and the Lehmann representations of correlation functions 
are ill-defined in the infinite volume limit as the normalization condition for scattering 
states is (for 9 1 > 9 2 > . . . > 9 n and 9[ > 9' 2 > . . . > 9' n ) 



O1 - an (0 1 , 



n \0' n , • • • , 9[) < .,. a{ = J] 2*5(9, - 9' 3 ) 5 aj>a ,. (12) 

3=1 

The idea of a low-temperature expansion is to subtract these divergences in some way 
[21, 23, 24]. Here we proceed as follows. We separate the contributions C r ^ s in the 
Lehmann representation of the two-point function (9) according to their different formal 
temperature dependencies into 

C£(«, q) = E°,(u, q) + q), (13) 

where 



n f d9\ . . . d9 r f d9\ . . . d9' 



P. 



X 



uo + iS - E s + E r 



a ^(9 1 ...9 r \O(0,0)\9' s ...9' 1 ) < .. < . (14) 
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The functions E r>8 , F r>s are related by 

e?M«) = [ F Z(-^-i)T- ( 15 ) 

The matrix elements ai - ar ($ 1 . . . 9 r \O(0, 0)\9' s . . . 9'i)a' s ...a', can be decomposed into a 
connected and a disconnected contribution. The latter is characterized by containing 
factors of 5(9 j — 9' k ), signalling that some of the particles do not encounter the operator 
O in the process described by the matrix element. A fundamental assumption of our 
approach is that the disconnected contributions act to cancel the partition function in 
the denominator of (9). More precisely we define quantities 

n 

rpO,n _ rpO _ \ ~~ 57 rpO,n-m 
111 j,k ~ ^ 'j,k 2.^1 771 j~m,k-m ' 
m=l 
n 

F?f = F° k -Y. Z rnF?^ m T m ,n = 0,1,2... (16) 

m=l 
n— 1 oo 

^ = E<^ + E S S ^ = 0,1,2... (17) 

k=0 m=n 
n—1 oo 

?n= » n = 0,1,2... (18) 

k=0 m=n 

The key assertion of our low-temperature expansion is that the quantities defined 
in this way are finite in the thermodynamic limit. Upon re-ordering of the infinite sums 
the two-point function (9) is expressed in terms of the £f and as 



xo\y> 



q) = Y;£?(u,,q) + F?{u,,q). (19) 



s=0 



By construction £f and Tf formally have a temperature dependence 

£?,T? ~0(e^ A ), (20) 

and (19) hence constitutes a low-temperature expansion of Xo{w,o)- The two-point 
functions we analyze in detail below have the symmetry 

Xo(u,q) = Xo(-^,-g), (21) 

which relates the positive and negative frequency regions. It is useful to combine £i and 
T{ into quantities that exhibit the same symmetry 

C?(uJ,q) = £?(uJ,q)+J?(u,q). (22) 

In terms of these quantities the spectral representation takes the form 

oo 

Xo{w,q) = Y,C?(u,q). (23) 

A key property of the Cf is that they are finite in the thermodynamic limit. We 
demonstrate this explicitly for the first nontrivial terms Cf and below and postulate 
that it is true in general. 
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2.1. Resummation 

Following the procedure set out above, the finite temperature Lehman representation 
of the particular two-point functions analyzed below can be re-expressed in the form 
(23), where the quantities C®(u,q) are finite in the infinite volume limit. However, in 
the cases we are interested in, the functions C®(u,q) are not uniformly small. In order 
to make this statement more precise let us denote the single-particle dispersion relation 
by 



e{q) = v/A 2 + v 2 q 2 . (24) 

We observe that as long as both uo ± e{q) ~ 0(1), the C® are of order C^e~^ rA j and 
hence (23) provides a good low-temperature expansion of the two point functions we 
are interested in far away from the mass shell. On the other hand, when we approach 
the mass shell we have 

Cf(c,g)oc(^-. 2 (g)) r+1 . (25) 

In order to obtain an expression for the susceptibility close to the mass shell we therefore 
need to sum up an infinite number of terms in (23). For the cases considered below, the 
zero-temperature two point function is of the form 

C o"V <?) = 7 rrj + • • • , (26) 

{uj + to) 1 — e z (q) 

where the corrections are negligible in the regime of temperatures and frequencies we 
consider. We then introduce a quantity E (oo, q) by defining 

Xo{uJ,q) = - r. 27 

The low-temperature expansion (23) for x ^^) then provides a way of determining 
low temperature approximations to TPiuo^q) in the following way. Assuming that a 
low-temperature expansion of the form 

£>,<?) = ££>,<?) (28) 

71=1 

exists, we may determine the leading term at low temperatures by expanding 

XoCw, q) = CoV, q) + [CJV, q)] V^, q) + [C?{u, q)f [tP(u>, q)]\... 

= C °(^, q) + [C°(u, q)] 2 Sf ( W , g) + . . . (29) 
and then comparing this expansion to (23). This gives 
Sf(^,g)=Cf(^,g)[C o (^g)]- 2 , 

= -C?(u,q) [Ef (u, q)] 2 + C 2 (uj, q) [C°(u : q)] ~ 2 , 
Z%(u,q)=... (30) 

We now turn to the implementation of the programme set out above to the case of the 
quantum Ising model. 
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The Hamiltonian of the transverse field Ising ferromagnet is given by 

H = J2-J« + i + ho- x n) (31) 

n 

where we take J,h > 0. The phase diagram of the model (31) is shown in Fig. 1. At 



°r>=o 



<<>±o 



h/J 



Figure 1. Phase diagram of the quantum Ising chain. At zero temperature 
ferromagnetic long range order occurs for h < J. 

zero temperature the quantum Ising model (31) exhibits two phases: 

(i) Ordered Phase: 

This phase occurs for h < J and is characterized by spontaneously broken Z2 
symmetry associated with long-range magnetic order along the ^-direction 

«> ^ 0. (32) 

(ii) Disordered Phase: 

This phase occurs for h > J. There is no spontaneous symmetry breaking and the 
only ordered moment is along the x-direction. 

The two phases are related by the Kramers- Wannier duality transformation [42], which 
provides a map of operators and therefore correlation functions [43, 44]. By virtue of 
this map it is sufficient to consider the regime h > J only. For the remainder of this 
paper we will concentrate on this parameter regime. 

3.1. Scaling Limit in the Disordered Phase 

The quantum Ising model (31) can be solved exactly [45, 46]. However, rather than 
analyzing the lattice model (31) directly we concentrate on the simpler scaling limit 
defined by 

J, h — > 00 , a — > 0, I J — h\ — A fixed , v = a \/jh fixed. (33) 
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Here a is the lattice spacing. In the scaling limit the lattice spin operators turn into 

1 /8 

continuum fields a z n — > Ca cr(x), while the Hamiltonian is expressed in terms of left 
and right-moving Major ana fermions 



%v — — — 

—{ipd x ip - ipd x ip) - iAipip 



(34) 



The excitation spectrum of (34) follows from a mode expansion of the Majorana field 
and is given by 



e(q) = VA 2 + v 2 q 2 - (35) 

In the disordered phase excitations can be thought of in terms of simple spin flips, as can 
be seen by considering the limit h ^> J. Let us denote the corresponding annihilation 
and creation operators by Z{9) and Z*(9) respectively. They fulfil the simple Faddeev- 
Zamolodchikov algebra 

A(9 1 )A(9 2 ) = SA(9 2 )A(9 1 ), 

A\9 1 )A\9 2 ) = SA\9 2 )A\9 X ), (36) 
A(9 l )A\9 2 ) = 2n5(9 1 - 9 2 ) + SA\9 2 )A{9 1 ), 

where the scattering matrix is S = — 1. The ground state is then defined as 

A(0)|O)=O, (37) 

and a basis of scattering states is given by 

\9 1 ,...,9 n )=A^9 1 )...A^9 n )\0). (38) 

Energy and momentum of the scattering states are by construction additive 

s 

E S (9 U ...,9 S ) = J] A cosh ^ , 
i=i 
s A 

P S (9 1 ,...,9 S ) = ^-sirih^ . (39) 
j=i 

In terms of the states (38) the resolution of the identity reads 

id= \O){O\+^J_^^^\0u---,0n){9n,...A\. (40) 

For the calculation of correlation functions the knowledge of the matrix elements or 
form factors of local operators is necessary. In the disordered phase the non- vanishing 
form factors of a z contain an odd number of particles and are given by [46, 47, 48, 49, 5] 

2n+l f) - f) 

(0| a(0, 0) \9 U ... , 9 2n+1 ) = i n a J] tanh -L—L. (41) 

i<j 

It is customary to choose the normalization of the field o~(x) such that 

lim (0\a(x)a(0)\0) = — — f, (42) 



\X 4 
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which implies that 

i 

, .4= 1.28242712910062... (43) 

We note that in this normalization the continuum field o~(x) is related to the lattice spin 
operator by crj — > 2 1 ^ 24 e 1 ^ s A~ 3 ^ 2 a^ 8 cr(x). In what follows we need more general matrix 
elements of the form 

(9[,...,9' k \a(0,0)\9 1 ,...,9 n ). (44) 

These can be calculated using crossing relations following [1]. The necessary identities 
are summarized in Appendix A. 



a 



2i2e 



1 A 3 



A 

v 



3.2. Spectral Representation of the Dynamical Susceptibility 

Our main interest is in calculating the retarded dynamical susceptibility at finite 
temperature, which is obtained by analytically continuing the Matsubara two-point 
function 

Xa(u,q)= f drdxe^ - iqx Xa(r,x) . , 



X a(r,x) = -(7>(r,:r)a(0,0)>. (45) 

In the basis of scattering states introduced above, the following formal spectral 
representation for the finite temperature dynamical susceptibility holds 

.. oo 

;fc( W , g ) = - . (46) 

r,s=0 

Here C r>s denotes the contribution with r particles in the thermal trace and s in the 
intermediate state and is of the form (10) without isotopic quantum numbers aj and a' k . 
As we have already stated, both the partition function and the Lehmann representations 
of correlation functions are ill-defined in the infinite volume limit by virtue of the 
normalization condition of states (12). Setting this issue aside for a moment, we may 
cast C° s in the form 

^vn- i \ f d9\ . . . d9 r f d9\ . . . d& . 



(2yr) r r! J (2tt) s s\ 

p —0E r p —f3E s 

-L——I——\{9 X . . . 9 r \a(0, 0)\9> . . . 9[)\\ (47) 

In order to implement the low temperature expansion we separate the C° s according to 
their (formal) temperature dependencies into E° s and Ff s following (13), (14). At zero 
temperature the spectral sum simplifies dramatically as only terms with r = or s = 
remain. 
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3.3. Zero Temperature Dynamical Susceptibility 

At T = the leading contributions to the dynamical susceptibility at low frequencies 



arc 



X*(u, q) « g) + F£ (w, g) + ^(w, g) + ^(w, g 

Here the 1-particle contributions are 



(48) 



va 



1 



e(q) u - e(q) + iO' 
va 2 1 



e(q) lu + e(q) + iO 
The 3-particle terms can be cast in the form 

E5 >3 (u,q)= [Ff t0 (-u,-q)Y 

* r d 9 ld 9 2 1 tanh 2 (^)tanh 2 (^)tanh 2 (^ 



(49) 



va 



A J 6(27r) 2 cosh^ 



iO — A ^7=1 cosn + ^0 



where 



#3 = arcsinh y— — sinh 6\ — sinh #2 J ■ 



(50) 



(51) 



We plot the real and imaginary parts of E ^(uj,q = 0) in Fig. 2. In order to plot these 
functions it is useful to separate off a dimensionful normalization factor 



iV n 



va 2 



A 2 



(52) 



We see that by virtue of the smallness of Eq 3 the dynamical response at low energies 
is dominated by the coherent single-particle contributions E% x + F[ . This yields the 





20 

to/A 



Figure 2. Real and imaginary parts of E^ 3 (uj 7 q = 0). 

following result for the temperature independent part of the spectral representation for 
the dynamical susceptibility 

2va 2 



C >,g)«£^,g)+F* (u;,g) 



{to + iO)* - e 2 (q) 



(53) 
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where e(q) = a/A 2 + v 2 q 2 . 



3.4- Infinite Volume Regularization 

At low temperatures the next most important contributions arise from Cf 2 and C 2 i> 
which are formally given by 

, dfldflA e- Af3c ^ - e -A^Kffi)+c(ft,)] 



2(2tt) 2 + ?0 - A[c(0i) + c(0 2 ) - c( 
x |(0|a(O)|0 2 , 9 1 )\ 2 S(vq - A[ S (^) + s(0 2 ) - s(0)]), (54) 

/• rfflrfflA e -A/3[c(gi)+c(e 2 )] _ e -A/3 C (0) 

2,i{ u ^)- v J 2 (2vr) 2 w + iO + A[ C (0 1 ) + c(0 2 )-c(0)] 

x |(0KO)|0 2 , 0i)| 2 <S(«9 + A[a(0O + s(0 2 ) - s(0)]), (55) 
where c(0) = cosh0 and s(9) = sinh0. We note that 6*2,1 can be obtained from C\ y2 as 
C^(u,q)= [CU-uj, -q)]*. (56) 

In order to proceed further we need to evaluate the products of form factors. Following 
Smirnov [1] we have the following crossing relations for three-particle form factors of 
the spin field in the disordered phase 

(0 1; 9 2 \a(0, O)|0 3 > = (0i - i0, 2 - z0|<r(0, O)|0 3 ) + 2tt^[5(0 32 ) - 5(9 31 )] 

= (0x + iO, 2 + z0|(t(0, O)|0 3 > - 2ira[5(9 32 ) - 6(9 31 )\ 

= (0x + zO, 02 - z0|<r(0, 0) |0 3 ) + 2na[5(9 32 ) + 5(9 31 )] 

= (0! - iO, 2 + i0\a(0, O)|0 3 ) - 27ra[<5(0 32 ) + 5(0 31 )], (57) 

where we have defined 9jk = 9j — 9^. The poles occurring in the form factors appearing 
on the RHS have been shifted away from the real rapidity axis. The delta-function 
contributions correspond to disconnected pieces of the form factor. It is clear from 
expression (57) that the absolute value squared of a form factor is ill-defined as a 
consequence of working in an infinite volume. As we will show below, the absolute value 
squared of a form factor contains divergent pieces that get canceled by a corresponding 
divergence in the partition function. In order to exhibit these cancellations we need to 
exhibit the divergences in the form factor squares explicitly. For the three-particle form 
factor of the spin field we do this as follows 

1(03^(0, 0)|0 2 ,0i)| 2 = lim(0 3 k(O,O)|0 2 , 9 1 )(9 1 , 9 2 \a(0, 0)|0 3 + «). (58) 

K— >0 

The product of form factors on the RHS is now well defined, so that we can use (57) 
to extract the divergent pieces explicitly. For form factors involving more than three 
particles one needs to introduce several parameters kj, e.g. 

|(0 4 ,0 5 |a(O,O)|0 3 ,02,0 1 )| 2 = lim (04,05^(0,0)103,02,0!) 

X (01,02,0 3 KO,O)|0 5 + «2,04 + Ki). ( 59 ) 
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(60) 



In order to calculate the contributions Cf 2 and C 21 by means of contour integral 
techniques the following choice of analytically continuing 12 into the complex plane 
is the most convenient 

\(9MoMWi)\ 2 = 

lim ({6 1 + zO, 9 2 - i0|<r(0, 0)|6» 3 + k) + 27ia[5(9 32 + re) + 5(9 31 + k)] 

x ((0|a(O, 0)|6> 2 - iO, a + iO) - 2trx[5(0 32 ) + 5(0 3 i)]) • 

The product of from factors on the RHS of eqn (60) is evaluated in Appendix B. Using 
(B.14) we can express the contribution C" 2 m the form 

Cf )2 (u,, q) = C$">, q) + Cg(w, g), (61) 

where the "connected" Cf^ 1111 ^, q) and "disconnected" Cf 2 (u;,g) parts correspond to 
the contributions of the first and second terms on the rhs of (B.14). The disconnected 
part is given by 

e -A/3c(e 3 ) _ e -A/3[c(0i)+ C (0 2 )] 



C™{u,q)=va 2 \\m / decide- 



W *0 - A[c(0 x ) + c(0 2 ) - c(0 3 )] 



X S( V q - A[ S (0!) + s(0 2 ) - S{B 3 )]) [S{k)S{B 31 ) + 5(032)5(031)] 



vex 2 1 



-/fe(<z) 



e(g) u; + iO - e(q) 



(62) 



We note that most of the "cross-terms" derived in Appendix B cancel in the integral 
as they are antisymmetric in 9\ and 2 , while the remaining part of the integrand is 
symmetric. The connected part of Ci |2 is given by 

d9d9+d9„ 



C c °™(uj,q)=va 2 



■A/3c(0) 



-2Ac(0+)c(0_) 



x K(9_,9+, 
where we have changed variables to 

a 9 2 ±e 1 

± 2 ' 
and defined a function 

K(9_,9 + ,9 



(2tt) 2 l ~ 1 
5(ug - A[2s(0+)c(0_) - s(9)]) 



uj- A[2c(0+)c(0_) -c 



tanh 



tanh^ 



7_-fl+fl + -i0 



tanh^ 



?_+e-e + -to 



(63) 



(64) 



(65) 



We can carry out the integration over 0+ using the momentum conservation delta- 
function. The result is 

d9d9^ e - A/3c ^ - e-ZMsAM K(6_, 9° + (q, 8, 0_), 0) 



C 1 T n (^,g) = w 2 



(2vr) 2 n(9,u)-u(q,e,e. 



M (g,0,0_ 



-.(66) 
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Here we have introduced the notations 

0,(6, uj) = uj + i5 + A cosh(6 l ) , 

Q(9,q) = g+-sinh(#), 

v 

"(9,9,9-) = V("Q(9,9)) 2 + (2Ac(0_))2. (67) 

In order to proceed it is useful to consider the integrand as a function of the complex 
variable (the analytic properties as a function of 6 are not as nice). We want to 
move the integration contour away from the real axis in order to avoid the vicinities of 
the double poles. The branch point of the square roots and inverse hyperbolic functions 
occur only at |Im(6L)| = |, which allows us to move the contour to a line parallel to 
the real axis in the lower half-plane. When doing this we may encounter a simple pole 
when 

n(e,w)-u(q,e,ej) = o. (68) 

Two kinds of solutions to this equation exist in the strip — | < 7 < Im(6L) < (for 
simplicity we assume uo > in the following) 

(i) If Q 2 — (vQ) 2 is larger than 4A 2 , we have a simple pole at 

a(uo,q,6) = — arccoshf^—^-^-") — iO , (69) 



2A 

s(u, q, 9) = [(lo + A cosh 9) 2 - (vq + A sinh 9) 2 ] * . (70) 
(ii) If 4A 2 cos 2 7 < Vt 2 - (vQ) 2 < 4A 2 , we have a pole at 



a(uj, q, 9) = —1 arccos 

Defining 



s{u,q,\ 

we may cast C^f 111 in the form 



@o(u,q,6) = arcsinhf ^ + AS 7n , (72) 



^ e -/3Ac(e) ^ 



CT° nI1 (u;,g) = -w^(l-e' 
- va 2 (l - e-P") 



+ 27rs(u),q,0)y/P(u},q,e)-4A* 
d9 e~ A ^ K(a, 9 , 9) 
T7 27r s(u, q, 9) v/4A 2 -s 2 (u,q,9) 



d9 f ^AcW-e-W^-ljf^,^,, 



+ W ./ (2Tr) 2 J s d9 -n(9,uj)-u(q,9,9-) u(q,6,6-) ' (?3) 
Here S is a straight line in the lower half-plane parallel to the real axis with imaginary 
part — 27, S + are the segments of the real axis such that s 2 > 4A 2 and are the 
segments of the real axis such that 4A 2 cos 2 7 < s 2 < 4A 2 . The segments S+, T+ for 
positive frequencies uj > can be characterized as follows: 
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e < 



(—00, 00) 

(—00, a_] U [a+, 00) 
(—00, a_] 
[a+, 00) 



if > |ug| and s 2 > A 2 
if uj > \vq\ and s 2 < A 2 
if < uj < v q , 
if < uj < —vq , 



(74) 



where 



(ii) T+: 



a± 



In 



3A 2 



s 2 ± Vs 4 - 10A 2 s 2 + 9A 4 



2A(cj - vq) 



(75) 



e < 



[a. 
[a- 



„/ 1 



U [a + ,a + \ 



[a' + , a+] 



if uj > \vq\ and A 2 7 2 < s 2 < A 2 
if uj > \vq\ and s 2 < A 2 7 2 , 
if < uj < vq , 
if < uj < —vq , 



(76) 



where j± = 1 ± 2 cos 7 (we have assumed that < 7 < |) and 



In 



-s 2 - A 2 7+7 _ ± y/(s 2 + A 2 7+7 „) 2 - 4A 2 g" 2 
2A(c^ - vq) 



(77) 



The real part of C^ ™ is given by the last two terms in (73). For fixed 7 there always 
will be a value of 9 such that the 6L integral in the third line of (73) is very close to a 
singularity at = 0. The problem occurs at 9 such that 

Q(9, uj) - V(vQ(9, q)) 2 + 4A 2 cos 2 ( 7 ) = 0. (78) 



For uj > \vq\ there are singularities at 9 



if s < A 71 and for uj < vq there is a 



single singularity at 9 
the regime 

s 2 « A 2 . 



From a practical point of view we are interested only in 



(79) 



For 7 = I we have 7_ = 0, so that we basically can always stay away from this problem. 
If we want to know the answer for s 2 < we can always calculate f(9,uj, q) for different 
values of 7 and in this way always stay away from having to deal with a singularity. 
This concludes our evaluation of 62,1(0;, q) in the infinite volume regularization. The 
contribution C 2> i(uj, q) can be obtained in the same way. 

We are now in a position to calculate the the leading contributions to the quantity 

C l ~ E l,0 + F 0,l + ( E l,2 ~ 2 l E 0,l) + ( F 2,l ~ Z l F lfi) ■ 

The leading corrections to Cf are (Ef 4 — Z\Eq 2 ) and (F^ r 1 — Z\F^ Q 
expected to be small for the same reasons that Eq 3 is negligible for uj s 
observe that F% x — Z\F° Q is negligibly small in the parameter regime we are interested 
in [uj A and low temperatures). We therefore drop it in the following. 



(80) 

and these are 
A. We further 
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As Z\ is ill-defined in the infinite- volume limit we regulate it by "shifting the trace" 
in the same way as (60) 



Zl = lim / — e-^(9\e + K) = \im5(K) [ d9 e~^ e \ 
k— >o J 2tt k— >o J 



51) 



Combining (22), (61), (62), (73) and (81) we arrive at the following expression for the 
first subleading part of the dynamical susceptibility (23) 



Ci(u, q) ~ — iva 



de e -^ Ac W K{a(u, q, 9), 9 {u, q, 9), 9) 
S+ 2^ s{u,q,9)y/W(^q,0)-l& 

d9 e~ A ^ e ) K(a(u, q, 9), 9 (uj, q, 9), 9) 
Tl 2^ s{u),q,Q)yJZE*-P{u,q,6) ' 

dQ , e -pAc(6) K (9_,9° + (q,9,9_) 
va I ~tz — rr: / gs6L 



v a 2 



(2tt) 2 7 s [n(e,u)-u(q,e,e-)]u( q ,e,e- 



(82) 



Here S + and are the segments of the real axis characterized by s 2 (u,q,9) > A A 2 
and 4 A 2 cos 2 7 < s 2 (u,q,9) < 4 A 2 respectively and S is the contour from —00 — 27 to 
00 — i'j parallel to the real axis. 

3.5. Finite Volume Regularization 

Another way to regulate infinities in matrix elements is to work in a large, finite volume 
R. The Hamiltonian on finite, periodic line of length R is 

R dx ~ 
2^ 

The Hilbert space of the theory divides itself into two sectors: Neveu-Schwarz (NS) 
and Ramond (R). The NS-sector consists of a Fock space built with even numbers of 
half-integer fermionic modes, i.e. states of the form 

\Pi"-P2n)ns = a p 1 -'- a p arf \ Q )NS , (84) 

where a mode's momentum satisfies 
2tt / 1^ 



H 







%v — — — 
—(ipd x ip - i)d x ip) - iAipip 



(83) 



Vi = ~R \ i + 2/ ' n ' 1 E Z ' ^ 
On the other hand, the R-sector consists of a Fock space composed of odd numbers of 
even integer fermionic modes, 

2ix 

\h--- k 2M +i)R = 4i • • • 4 2M+1 » k i = -R n i- (86) 
Energy E(pi) and momentum P(pi) of a NS state \pi ■ ■ -P2n)ns are given simply by 

2N 2N 

E({ Pi }) = j2 £ &)> p (fe» = E^> ( 8T ) 

i=l i=l 
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where as before e(p) = y A 2 + v 2 p 2 . An identical relation holds for states in the R 
sector. It is useful to parametrize the momenta in terms of a rapidity variable as 
A 



Pi 



sinh(0 Pi ) . 



The finite volume form factors of the spin field have been determined in Refs [26, 50]. 
The non-vanishing form factors in the disordered phase are 

,-M+AT. 



R(h---k 2 M+i\o-(0)\pi-- 

' Oki — Ok, 



J tanh 



P2N)NS = I 

J tanh 



n , — 



'■J 



Pj 



Jcoth 



Ok, — 0. 



(89) 



where for large AR we have 

C R =a + o(e- AR 

1 



i<j 



i -J 



O e 



9(0) 



AR 



(90) 



'ARv~ l cosh(fl) 

We note that the factors g(0) disappear when the states are normalized in terms of 
rapidity variables. Importantly, up to exponentially small corrections, the matrix 
elements (89) have the same functional form as at R = oo. The essential difference 
between large but finite R and the infinite volume limit is that at finite R, the momenta 
are quantized. This is a pattern that repeats itself for general integrable models, as 
emphasized in Ref. [25], and that we will exploit for our analysis of spin-1 chains. 

Crucially, for large but finite R the matrix elements (89) are finite. In the infinite 
volume limit divergences develop in the factor 



Jcoth 



Ok, - o Pj 



(91) 



1,3 



and occur when two momenta, ki and pj, approach one another. However, a finite R 
regulates these divergences by virtue of ki lying in the R-sector with integer quantization 
and pj in the NS-sector with half-integer quantization respectively. Hence the two are 
never exactly equal. The finite temperature Lehmann representation for the two-point 
function of the spin field on a ring of length R takes the form 



1 oo 



u,q) 



(92) 



r,s=0 



where e.g. 



2-nn q 
~~R~ 



R 



dx 



aire 



iu) n r—iqx r~iR 



{ k j}Api\ 

e -f3E({k J }) _ e -(3E({ Pt }) 

X u + irf-EdpiV + Edkj}) 



■P2n) 



TVS I 



|2 



R5 



P{{Pi})-P{{kj}),q > 
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oo 

Z R = \ + Y, e ~ Pe{p) + E e~^ +£ ^ + ... = J2 Z n- ( 93 ) 

pG-R. pi,p 2 GNS n=0 

Here 77 is a positive infinitesimal. All terms in the expansion (93) are finite. As before, 
we re-order the spectral sums according to (22), which gives 

00 

x^<?) = E c ^)> ( 94 ) 

r=0 

where the C R are defined as the finite volume analogs of (22). 

3. 6. Comparison Between Infinite and Finite Volume Regularizations 

In this section we establish the equivalence between the infinite volume regularization 
used in Section 3.2 to evaluate the leading order terms of Xafo, q) and the finite volume 
regularization (as R — > 00) introduced in the preceding section. In particular we 
establish that C R + C R , the first two terms in the temperature expansion of x R { u) i ( l)i 
and given by 

C R s=! F, R 4- F R 

C R « E R + F R + {E R 2 - Z R E R :i ) + {F R i ~ Z?F*) , (95) 

are equal to their counterparts in infinite volume once we take volume R to infinity. 
The equivalence of Co in both schemes is straightforward. We have for C R ± and C R the 
following 



C^ir, x) = £ «(t, x) + F *(t, x) = -a 2 v ^ 



e -e(pi)r+ipx 



p -£{Pi)(J3-r)-ipx 

C R (r, x) = E r (t, x) + F R (r, x) = -a 2 v £ . (96) 

p% R£ ^) 



Taking the Fourier transform we find 



+ tO - e(q) ' u,iV ' e(q) lu + iO - e(q) ' 

p , a 2 w e-^W ^p , . a 2 v -1 

E R (LU,q) = —— — , Ffn(u,q) = — — 7-r. 97 

1M e(g)^ + z0 + e(g) > ^ ^ + z0 + e(g) v ; 

In comparing this to (49) and using F R s (u>,q) = [E R r (— u>, —q)]* we see that the finite 
and infinite volume regularizations for Cq lead to the same result. 

The above also establishes that the first two terms of Cf , E° + Fq ± , in the two 
schemes are the same. It leaves then to verify that 

Mm (E R 2 - Z R E R A ) = El 2 - Z 1 ES tV 

lim (F R - Z R F R ) = F^ - Z X F^. (98) 
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We first consider (E^ 2 ~ ^l'^oi) ■ m a finite volume this term is equal to 

rrR r~*R. 0~ 2 V f ^ — > ^ — > 



^ ^ ; C ^ ^ ^ ^ ^g,p2+P3-pi 

p l6 * P2>P3 eNS UJ + l0 ~ £ ^ ~ £ ^ + £ M 
e -/fe(pi) tanh 2 (^f) coth 2 (^) coth 2 (^f) 

e(Pi)e(P2)e(Ps) 
_ ^ y ^pl y e -*(*0 (99) 



We rewrite i^f* using the identity 

Z n = J^^W^^l^ 



7T 2 ^ (71+ 1/2) 2 

E E jpfo-p,*. - < 100 > 



Inserting this into eqn (99) we obtain 

Pl&R P2,P3&NS L 



} q,P2+P3~Pl 



u + iO- e(p 2 ) - e(p 3 ) + e{p l ) 



e -/fe(w) tanh 2 (^f) coth 2 (^f) coth 2 (^) 

e(Pi)e(P2)e(ps) 

_ 1 <W% 4erf^ 

e(Pa + I)) w + iO - e(p 3 + |) v 2 ( Pl - p 2 f 



1 <W+- 4e-/3 £ (pi) 



e (pa + I)) oo + iO - e(p 2 + f ) w 2 (pi - p 3 ) 2 



(101) 



We now see that the above expression has no singular pieces of the form (pi — p 2 )~ 2 or 
(pi ~ Vz) 2 as Pi approaches p 2 or p 3 . There are singular terms of the form (pi — p 2 )~ l 
and (pi — P3) _1 , but these can be seen to vanish upon summation and so can be ignored. 
Given the absence of singular terms we can take the limit R — > oo and convert the 
summations to principal value integrals, i.e. 



— 2 C 
87T J 



5( vq -A(s(e 2 ) + s(6 3 ) -s m) 0J3 cQth2 ^ cQth2 ^3 

u + iO - e{6 2 ) - e{6 3 ) + e^) y 2 ! y 2 ! y 2 1 

5(vq - As(6 3 )) te^efa) 5(vq - As(6 2 )) 4e(0 1 )e(0 3 ) 



u + iO - e(6 3 ) A 2 (s(9 1 ) - s(6 2 )) 2 u + iO - e{6 2 ) A 2 ( S (^) - s{6 3 ))\ ' 
The final step before being able to compare to the results stemming from our infinite 
volume regularization scheme is to convert the principal value integrals into integrals 
along contours deformed by iir/'s about the singularities found at 9\ = 9 2 and 6\ = 9 3 . 
To do so we employ the identities 



) ' (3(6,) - s(6 2 )y J Cv{d2) 1 (aft) - 



Finite Temperature Dynamical Correlations in Massive Integrable QFTs 



22 



d9 x 



m- 



8\e 



21 j 



2 5(9 



12, 



and 



f (01,02 



M2 tanh 2 ^ 



old. 



c{v 2 )c(0i) vc 2 (e 2 ) 

I (01,02) 



f(0i,0 2 ) + O(r]), (103) 



C v ( 6l ) 2 tanh 2 (^) 



d0 2 f (0i,02) 



Am5'(6 2l ) + -8(9 l2 ) 
V 



+ o( v ), 



(104) 



where f(6i,6 2 ) is a test function and the integration contours are defined in Fig. 3. 
Deforming the contours into the lower half plane instead changes the signs of the 8'- 
terms. Using (103) and (104) to deform the #3 and 9 2 integrals into the upper and lower 



Im 9? 



Im 9 P 



Re9 P 



Re 9? 



e r Ti a, -hi 



C n (0,) 



Figure 3. Integration contours R^O) and C v (6). 

half-plane respectively, we can rewrite E R 2 — Z R E R X in the following way 

E R - Z R E R 
^1,2 ^1 -^0,1 



va 2 



SIX' 



tanh 







23 



coth 2 f^^)coth 2 (^^ 



ou + iO-e(9 2 )-e(9 3 ) + e(9 1 ] 
) + 16n 2 5'(9 12 )5'(9 13 ) 



+ 4mC oth 2 ( 9 -^±^l)5'(9 13 ] 



A7ricoth ^^_jR)S'(9 12 ) 



va 2 
2^2 
va 2 



- / 3e(6 1 ) S ( V( l ~ As ( d 2)) 



6(01)6(03 



d9 1 d9 2 d9 3 e Lu + iQ _ ^ _ _ 

d9 ld 9 2 d9 3 e-^ 6 ^T As ^) 



2tt 2 J u + i0-e(9 3 ) A 2 (s(9 1 + ir ] ) - s(9 2 )) 2 ' 

We note that all terms singular in rj that come about from deforming the contour of 
integration of 9\ about the singularities at 9 2 and 9 3 vanish (as they must as the principal 
value integral is well defined and finite). The last two integrals in (105) vanish as may 
be seen by deforming #3 — > 9 3 + in/ 2 and 9 2 — > 9 2 — in/ 2 respectively. After carrying 
out the integrals over the derivatives of delta functions we arrive at 



Cj \,i ^1 -^0,1 



Rt?R 



va 



w+i0-efl 2 )-e(»3 
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x 



tanh 2 (^) coth 2 (^p) coth 2 (^±^) + 8n*6(9 12 )6(9 13 ) 



(106) 



We are now in a position to show that E^ 2 — Z^Eq X equals its value in our infinite 
volume regularization scheme. Indeed, combining (63), (62) and (81) and then keeping 
only the part of order C(e _/3A ) we obtain (106). 

By an analogous consideration we can show that lim^oo \ — Z^F^ Q recovers the 
result of the infinite volume regularization scheme. 

3.7. Resummation 

Above we have argued that the finite temperature Lehmann representation can be 
reexpressed as 

oo 

Xa {uj,q) = Y,C:{uj,q), (107) 

where the quantities C^(u, q) are finite in the infinite volume limit. In particular 

C o 9) = ( T 2? n + " " " ( 108 ) 

(u + id) 2 — e 2 (q) 

is the zero-temperature dynamical structure factor. We emphasize that we only 
consider frequencies close to A and low temperatures, so that we can neglect n-particle 
contributions to Cq with n > 3 as their contributions are vanishingly small. For 
uj — e{q) ~ 0(1) C° is of order 0^e~ l3rA j and hence (107) provides a good low- 
temperature expansion of the dynamical susceptibility far away from the mass shell. 
On the other hand, when we approach the mass shell we have 

C s >,g)oc (u 2 -e 2 ( q )Y +1 . (109) 

In order to obtain an expression for the susceptibility close to the mass shell we therefore 
have to carry out a resummation of terms. This is achieved by introducing a quantity 
E(a>, q) by (27) and using our results for Cf and C% to carry out a low-temperature 
expansion of H(uj,q) using (30). 

4. Results for the Low- Temperature Dynamical Susceptibility of the 
Quantum Ising Model 

In order to present results obtained from the resummation (27) it is useful to define 
quantities 

>)( w n) = C o( u >l) (UQ) 

Xa [ ,9) l-Cg(u,q)E7 =1 ^)(u;,qy liiUj 

where Ti^\u,q) are given by (30). Loosely speaking y^a\u,q) is the dynamical 
susceptibility obtained by calculating the H{u,q) up to order C^e~( ra+1 ^ A J . From 
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the point of view of applications to experiment the relevant quantity is the dynamical 
structure factor, which is related to the retarded susceptibility by 

S(u, q) = — - — ^— ^Im xA u > <?)• (HI) 
7r 1 — e t 

In analogy with the definition of the m th order low temperature approximation x^ for 
the dynamical susceptibility we define 

S^(u,q) = -1— L^Im x^^.q)- (112) 
7r 1 — e t 

The leading order result S^ 1 ' (u, q) is most easily calculated using the infinite- volume 
regularization scheme and carrying out the integrals in (82) numerically. In order to 
determine S^(cu, q) we employ the finite volume regularization scheme instead. We 
calculate S™(u + irj, q) for system sizes up to R = 800. The non-zero imaginary part is 
introduced in order to suppress finite-size effects. If 77 is sufficiently large the numerically 
accessible values of R will coincide (within numerical accuracy) with the thermodynamic 
limit results. In order to obtain results for real frequencies we compute S^ 2 '(uj + irj, q) for 
a sequence of different 77's and then extrapolate to 77 — > 0. We have tested this method 
for S^'(u, q), which we can calculate directly in the thermodynamic limit, and found it 
to work well. 

At very low temperatures it is sufficient to calculate the leading approximation 
(lu, q) as the corrections are exponentially small in e~ A ^ T . We find that the difference 
between S^(u,q) and S^(cu,q) is negligible for T < 0.3 A and we discuss this regime 
first. Evaluation of S^'(uj, q) for low temperatures shows that, as expected, the T = 
delta function at u> = e(q) broadens with temperature. We find that the resulting peak 
scales as 

A /A\ 
peak height oc — exp^-j, 

peak width oc — expf— -Y (113) 



A V T 

In order to exhibit the evolution of the structure factor as a function of frequency for 
fixed momentum q with temperature it is therefore useful to rescale both the frequency 
axis and the structure factor. The result is shown in Fig.4 for a temperature range 
0.15 A < T < 0.3A, which corresponds to approximately a factor of 50 difference in 
peak height (and width). We see that the lineshape is asymmetric in frequency, with 
more spectral weight appearing at higher frequencies. The asymmetry increases with 
temperature. This effect is most easily quantified by comparison with a Lorentzian 
lineshape, which is done below. 

If we increase the temperature beyond T « 0.3A the correction T,^(uj,q) (30) 
in the expression for the dynamical susceptibility (27), (30) is no longer negligible 
compared to the leading term SW(o;,g). In Fig. 5 we show the leading S^(u,q) and 
improved S^(u,q) low-temperature approximation to the dynamical structure factor 
for T = 0.45A and T = 0.5A. As expected the difference between the two is small 
sufficiently far away from the mass shell. The most important effect is the shift of the 
maximum to higher frequencies. This effect is known as "temperature dependent gap" . 
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- T=0.25A 




T=0.3A 




(A/T) e A/T (co/A-l) 



Figure 4. Leading low-T approximation (id, q = 0) of the dynamical structure 
factor for several temperatures. Both axes have been rescaled as discussed in the text. 




3 4 

55 



T=0.5A 




1 

to/A 



Figure 5. Leading S^(uj,q = 0) and improved (u>, q = 0) low-T approximation of 
the dynamical structure factor for temperatures T = 0.45A and T = 0.5A respectively. 
The most important effect of the next to leading order approximation is the shift of 
the maximum upwards in frequency ( "temperature dependent gap" ) . 



4-1. Comparison to Semiclassical Results of Sachdev and Young 

In [28, 38] Sachdev and Young have developed a semiclassical approach to determine 
the dynamical structure factor. Their result §, valid for T A, is 

dt 



2n 



dx e iujt - iqx K(x,t)R(x,t) 



(114) 



Our definition of the structure factor differs from [38] by a factor of 2ir. 
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where 



a 2 



K( x ,t) = -K (\AW(x/cf-t 2 ) , 

7T 

R{x,t) = exp (- J°° ^ e - £ ^ T \x - v(k)t\ 
_de{k) v 2 k 

v ( k ) = —nr = 77TX- ( 115 ) 



Here a is defined in (43). At sufficiently low temperatures and when \vq\ <C y/TA the 
semiclassical result is well approximated by a Lorentzian [38] 



a 2 v 



S Loi (uj,q) = — --- (116) 

7T £ (g) (u - e{q)) 2 + 1/rj 



where 



r* = ^ . (117) 

We expect our low-temperature expansion to reproduce the semiclassical results at 
sufficiently low temperatures T<A. The comparison of the improved result S^(u, q) 
to the leading order S^(uj, q) suggests that the latter is a good approximation up to 
temperatures of T m 0.3A. In Fig.6 we show a comparison of our low temperature 
approximation to the semiclassical result (114) and the Lorentzian approximation (116) 
for T = 0.3A. We see that while the gross structure of the lineshapes is quite similar, 
there are considerable differences away from the mass shell. This implies that the 
semiclassical approximation is no longer quantitatively accurate for T = 0.3A. 




Figure 6. Dynamical structure factor for T = 0.3A, q = 0. The result of the present 
work (red line) is still in good agreement with the semiclassical result of Sachdev and 
Young as well as the Lorentzian approximation to the latter. The asymmetry of the 
lincshape is more pronounced than the one predicted by the semiclassical result. 

In order to exhibit the differences between the various approximations more clearly, 
we plot the ratios S^(u, q) / Sl oi (uj , q) and S sc (u, q)/S^ OT (uj, q) in Fig.7(a) for T = 0.2A. 
We see that the semiclassical approximation underestimates the asymmetry of the 
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lineshape. On the other hand, for sufficiently small temperatures the low-temperature 
expansion indeed recovers the semiclassical result as can be seen from Fig. 7(b). 




Figure 7. (a) Ratios of the dynamical structure factors calculated from eq. (27) 
(this work) and (114) (Sachdcv and Young) to the Lorentzian approximation (116) 
for T = 0.2A and q — 0. The non-monotonic behaviour very close to the mass shell 
indicates a small shift of the maximum of S^'(u), q — 0) away from the T = gap. (b) 
Ratio of SW(w,q = 0) to the Lorentzian approximation (116) for T = 0.1A and q = 0. 
We see that the low temperature approximation recovers the semiclassical result at 
sufficiently low temperatures. 

The ratio S^'(uj, q)/Si jOT (uj, q) displays non-monotonic behaviour close to the mass 
shell. This suggests that the maximum of S^'(uj, q) is in fact very slightly shifted away 
from the T = gap. This is indeed the case. 



4-2. Comparison to Exact Finite Temperature Form Factors 

In Refs [24, 27] the following form factor expansion for the two point function of the 
spin field in the quantum Ising chain was derived 

xAt,x) = yj£>»,(T,x) , (118) 



^.(r,x)=-^e-(«M/n 

J 7_1 



dBj 
2^ 



_f(0 ) e rM*i)-i#"( 9 J>-»H-(«i) 



)Ac(9' l )+i^s(e' l )x-r,-(9 l 



I) f[ tanh 2 (^p 



n>m=l 



X 



^ tanh 2 (^)ncoth 2 ^" "^ + Z ° 



p>q=l p,n 

where c(6) = cosh^, s(9) = sinh^, = 9j — 9^ and 



(119) 



C(P) = crexp 



(A/5) 2 [d9 l( 19 2 s(9 l )s( 



2 J ( 27r ) 2 s(a/3c(9 1 ))s(a/3c(9 2 ) 



In 



coth 



'12 
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KHJ irv ^ 2k - 

k=l 

ooq=iO 



1 



Ki{^2k - 1)PA 



±2 



ooq=iO 



V±(0) 

At low temperatures we have 
C((3) « a, 

e v±(e)~ [1 - 2e^ Ac ^] exp 



dff 



2vri sinh(# - & 



2% 
T— 

7T 




dx 



sinh x 



-PAc(0-x) 



-f}Ac{9+x)\ 



n(P) 



-/3A 



2A 

7ij3v 2 



An important question is in which way this expansion is related to the low-temperature 
expansion based on the zero temperature form factors. Fourier transforming and 
analytically continuing to real frequencies we find 

dO e-"-W 2n{(3) 



DUu;, q ) = C 2 ({3) 



DUu,q)= -C\(3) 



2ix u + iO - Ac(0) n 2 (/3) 



2n(/3) 



2tt w + «0 + Ac(0) n 2 (/3) + (g + fs(0)) 2 ' 

The main difference between D r s and the corresponding zero temperature quantities is 
the replacement of the momentum conservation delta function by a Lorentzian of width 
n(f3). We note that the imaginary part of Do,i has a square root divergence for uo — ► A 
and it is necessary to sum an infinite number of terms in (118) to get a meaningful 
answer [24]. In order to recover our low-temperature expansion from the Lehmann 
representation in terms of exact finite temperature form factors we should expand the 
Lorentzian expressing approximate momentum conservation under the integral, e.g. 



2n((3) 



q\9) 



where q{6) 



n 2 (f3)+q 2 (6) 



2\2 



2*6{q{0)) + 4n(/3)- {q 
-Ann 2 (P)5"(q(9)) + 0(n 3 (P)) 



5. 0(3) Nonlinear a-Model 

We now apply the methods outlined above to the 0(3) nonlinear sigma model. Unlike 
the quantum Ising model the sigma model describes a strongly interacting theory 
featuring dynamical mass generation. The Lagrangian of the sigma model is given 
by 



C 



1 



dx 



-d t n ■ d t n - vd x ii ■ d x ii 



(121) 
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The 0(3) nonlinear sigma model describes the scaling limit of integer spin-S 
Heisenberg models [51] 

H = J y £ t S j -S j+ i, S 2 j = S{S+1). (122) 

j 

The velocity and coupling constant of the sigma model are related to the lattice model 
parameters by 

9 = | , v = 2JSa , (123) 

where ao is the lattice spacing [52]. The lattice spin operators, Sj, are related to the 
continuum fields by 



Sj ~ S{-l) j n{x) + l(x) , x = ja , (124) 

vg v / at 



where l(x) oc — n(x) x ^7^. The dynamical susceptibilities of the lattice model are 



given by 

xt^.Q) = ~ /^^E ^(T r ^ +1 (r)^) . (125) 

Substituting (124) into (125) we see that in the vicinity of the antiferromagnetic wave 
number (Q = ^- + q with \q\ <C the lattice susceptibility at low energies uj <C t'/^o 
can be expressed in terms of the two-point function of the n-field 

Xfat(^, Q) oc - / dr / e^" T "^(r T n z (r,x)n 2 (0)) . (126) 



kJn — ><5 — IUJ 

The n-field needs to be renormalized and in an appropriate scheme is related to the spin 
field <fr by 

n a = C$ a . (127) 
In what follows we analyze the two-point function of the sigma model spin field 

*>,?)= - r^xe^-^(T T $ a (r,x)$ a (0,0)) , (128) 

JO OJn^S — lLU 

keeping in mind that at low energies we have 

tf* t L,- + q) oc X >,g). (129) 

The 0(3) nonlinear sigma model is integrable [53, 54, 55] and the exact spectrum 
and scattering matrix have been known for a long time. The elementary excitations of 
the sigma model are a triplet of massive particles with spins S z = ±1, 0. It is useful to 
parametrize their energy and momentum in terms of a rapidity variable 9 

e(9) — A cosh 9 , p(6) = — sinh 9. (130) 

v 

It is convenient to choose a basis such that the spin operators S a act on single-particle 
states as 

S a \9) b = ie abc \9) c , a = 1,2,3. (131) 
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A convenient basis for the Hilbert space is formed by scattering states of these 
elementary excitations. In order to describe it one introduces creation and annihilation 
operators Z a (6) and Z\{9) respectively, fulfilling the Faddeev-Zamolodchikov algebra 
(4). Here a, b are 0(3) quantum numbers and S is the exact two-particle scattering 
matrix [54] 

S c at(0) = (Ji{e)5 ab 8 cd + a 2 (6)5 ac 5 M + cr 3 (6)5 ad 5 bc , 
2niQ 



{d + m)(e-2m) 
— in) 



+ m)(e-2ni) ' 



2iri(0 — in) 

= " (o + ,,)(o-Ly < 132 > 

We note that the S-matrix is a solution to the Yang-Baxter equation and fulfils 

[SSWY = S5(-&). (133) 

Using the Faddeev-Zamolodchikov operators, a Fock space of states can be constructed 
as in (5), (6). Energy and momentum of are by construction additive and given by 
(7). In the basis of scattering states introduced above, formally the following spectral 
representation for the finite temperature dynamical susceptibility holds 

.. CO 

r,s=0 

Here C r>s are given by (10) 

C* s (cu,q)= I dr I dxe iuJnT - iqx C* s (r, 

(«■ <A — f — " " " — f d0 l - ■■ -(3E r -T(E s -E r )-i(P r -P s )x 

r ' A ' ] - J (27r)T! J (2irYs\ 

x E E kMOi-.-e r moM^--A)v m ^- (135) 

bi,...,b r b' 1 ,...,b' s 

The partition function can formally be expressed as in (11). The low-temperature 
expansion is constructed by following the steps set out in section 2, e.g. the functions 
Ef s used in the decomposition (13) of Cf s are 

. . v — -\ v — -\ f d9-i . . . d9 r f d0\ . . . d& . „ , „ „ 

0*»> = E E / / -^r^i + p < - p > 

b U ...,b r &'..,&' J 



X 



e 

x 



l3E r 

F - F kMOi ■ ■ ■ v r mo, o)\e' s . . . e[) b ,... b[ \ 2 , (136) 

uj + to — h s + h r 1 

where E r and P r are defined in (7). Our aim is to determine the functions Ef , Tf and Cf 
defined by (17), (18) and (22) respectively and to use them to obtain a low-temperature 
expansion (23) of the susceptibility. 
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5.1. Zero Temperature Dynamical Response 

The low energy dynamical response in the 0(3) nonlinear sigma model has been worked 
out in [56, 57]. The zero-temperature form factors in the infinite volume have been 
determined by several groups [6, 7, 8]. The one and three-particle form factors of the 
spin field are [7] 

m a (O)\0)b = S ab , (137) 

7T 3 



<0|$ Q (0)|£ 3 A A>a 3 a 2ai = yV>(W(W(02l) 

^aai5a 2 a 3 ^23 + ^aa 2 ^a 1 a 3 (#31 — 2lTl) + S aa3 S aia2 Ri2 



:i38i 



where 



ft — ?7T R 

*W = 5pS^ taBh V (139) 

Note that these differ slightly from [7] because we use the a different normalization 
condition for the scattering states 

a (R\R') b = 27r5 ab 5(R-R'). (140) 

Some useful identities involving the function i/)(R) are 

< / n . \ ,,. m #coth 2 (#/2) 

ft A + 7T Z 

^(ft + m+) =i)(R + m~) + — 6(9). (141) 

7T 

Here ir ± = n ± 0. 

At T = the leading contributions to the dynamical susceptibility at low 
frequencies are 

*>, g) » [J5*(o,, ?) + *J (f , <?) + JSjs(«. 9) + *?o(w, ?)] • ( 142 ) 

T=0 

Here the 1-particle contributions are 

= [*-(-*, -,)]* = ^ „,_^ ) + <Q . (143) 
The 3-particle terms can be cast in the form [57] 

2w 6 f°°dydz f(y + z)f(y-z)f(2z) [3vr 2 + 3^ 2 + y 2 ] 



[™ dy dz fjy + z)f(y - z)f(2z) [3tH + 3^ + yj 
7.^ (2vr) 2 s 2 - A 2 - 4A 2 cosh(^) [cosh(^) + cosh(y)] ' 1 } 



where 



/( g ) = ^+^) [tanh(z/2)]A ' 
s 2 = (w + iO) 2 - w 2 g 2 . (145) 
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We plot the real and imaginary parts of E® 3 (uj,q = 0) + F 03 (u),q = 0) in Fig.8. We 
see that by virtue of the smallness of Eq 3 + Ff the dynamical response at low energies 
is dominated by the coherent single-particle contributions Eq X + F* . This yields the 
following result for the temperature independent part of the spectral representation for 
the dynamical susceptibility 

C?(u, g) « Eg^v, q) + i£ ( W , q) = j- - ^ _ , (146) 



where e(q) = a/A 2 + v 2 q 2 . 

5.2. Infinite volume regularization 

At low temperatures the next most important contributions arise from Ef 2 and F 2 \, 
which are formally given by 

dede 1( ie 2 e -P*m 



2(2vr) 2 u + iO - A[c(0i) + c(0 2 ) - 
x |6(0|$ a (O)|0 2 ^i) febl | 2 5N- A[ S (^) + S (^)- S (0)]),(147) 



if 2,l( W »?) = -V 



6 ^ 2 7 2(2tt) 2 a; + iO + A[c(^) + c(0 3 ) - c(0)] 

x | 6l62 (0i, e 2 |$ a (0)|e) 6 | 2 5(^ + A[ S (^) + s (0 2 ) - s(0)]), (148) 

where c(0) = cosh0 and s(8) = sinh0. In order to proceed further we need to 
evaluate the absolute value squares of form factors. The individual form factors can 
be analytically continued following Smirnov [1], see Appendix A for a summary. The 
various ways of analytically continuing the 3-particle form factors are 

blb2 (9x,9 2 \<s> a (o)\9 3 ) b = bib2 (e 1 + to 1 e 2 -io\^ a (o)\e 3 } b 

+ 2n5(9 32 )5 bb2 5 abl + 2ir8(9 3 i)6 a b 2 6 bbl 
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= blba (e 1 -io,e 2 + io\$ a {o)\e 3 ) i 

+ 27rf(fo)5& a (0 ia ) + 2^(^31)^(^12) 
= 6 lb2 {^l + z0,6> 2 + z0|^ a (0)|^ 3 >6 

+ 27^(032)^(012) + 2ir5(6 31 )5 aba 5 Vbl 

= 6 1 6 9 <0i-iO,0 2 -iO|$ a (O)|0 3 ) ll 

+ 2n5(6 32 )5 abl 6 bb2 + 2^(^31)^(012). (149) 

Clearly the absolute value squared of (149) is ill-defined as it contains squares of 
delta functions. We can circumvent this problem by introducing infinitesimal shifts 
of rapidities in one of the form factors. In the case at hand this leads to 

| M2 (0 1 ,0 2 |$ a (O)|0 3 ) b | 2 = lim 6(031^(0)102,0!)^ M, 2 (01,02|<nO)|0 3 + K;>6 

6l62 (01 " *O,0 2 + *O|$ Q (O)|0 3 > 6 * + 2^(032)^(021) + 2^(^1)5^(021 

blb2 (0! + iO, 2 - iO|$ a (O)|0 3 + n) b + 27r5(^ 32 + K)5 bb2 5 abl + 2ir5(6 3l + K)5 ab2 5 bbl 

(150) 

where in the second line we have used (149). Using the results of Appendix C we can 
separate (150) into a connected contribution r conn given by (C.3), (C.7) and disconnected 
contributions r dls>1 and r dls ' 2 given by (C.4) and (C.5) respectively. Concomitantly E\ i2 
can be split into a connected and a disconnected part 

^(w, q) = ^iTV' q) + E^(u, q). (151) 
The connected part is 

d0id0 2 d0 3 e -/3Ac(0 3 ) 



lim 

K->0 



2(2tt) 2 u + iO- A[c(0i) + c(0 2 ) - c(0 3 )] 
x 5(vq - A[s(0!) + s(0 2 ) - s(0 3 )]) L(0x, 2 , 3 ), (152) 

where 



("is + i0)l^v4^i3 + »<h (9 23 - ?0) 2 „ tv4 ^23 - iO 

2)2 L 

The disconnected term, E^f is equal to 



x (^-^2 ™th (^^— j coth ( — _ ) . (153) 



dd x dd 2 dB 3 e -/3A C (e 3 ) 



2(2tt) 2 lu + iO- A[c(0 x ) +c(0 2 ) -c(0 3 )] 

2 

x*(ug-A[s(0i) + s(0 2 )-s(0)]) ^if^AA), (154) 

i=i 

where r d _ ls '-' are given by (C.9) and (C.ll) respectively. Upon evaluation, this term 
simplifies to 

5(vq - A[ S (0i) + s(0 2 ) - s(0 3 )])e^ Ac ^) 



^! di 2 S (^,g) =u y d9 1 d9 2 d9 3 



u + iO- A[c(0i) + c(0 2 ) - c(0 3 )] 
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(02 2 + 7r 2)(02 2 + 47r2) 



+ 



where Z\ is defined below in (161). Note that unlike the Ising model, the finite 
disconnected terms here involve derivatives of 5-functions. 

To evaluate the connected term, El° 2 nn , we now proceed in complete analogy with 
the Ising case. We change variables to 9± = {9 2 ± 6' 1 )/2, carry out the # + -integral using 
the momentum conservation delta-function and then shift the integration contour 
down in the complex plane. This results in 

old e^ Ac W K SM (a(cu, q, 9), 8 (u, q, 9), 9) 



E™P(u,q) 



IV 



2tt 



S\UJ 



s 2 (w,g,i 



4A 2 



d9 e~ A ^ K SM (a(u, q, 9), (u, q, 9), 
Tl 2tt s(uj, q, 9) ^/4A 2 - s 2 {u, q, 9) 



+ v 



al9 



dd. 



e~^ K SM (9-,9° + (q,9,9„),< 



[Q(9, uj) - u(q, 9, 9-)} u(q, 9, 9.) ' ^ 

where £1(6, uj), d+(q, 9, 6L) and u(q, 9, 9 J) are given by (67), a by (69), s(uj, q, 9) by (70), 
a by (71) and 

K SM (9 + , 0_, 6 3 ) = L(9 + -d_,6 + + 6_, 9 3 ). (157) 

The remaining integrals in (156) are easily evaluated numerically. 

We can further simplify the disconnected terms by performing some of the integrals. 
Doing so we obtain 

v 35(k) 



K>A) 



e(q) uj + iO - e(q) 

v 1 
e(q) uj + iO - e(q) 

v A 
e(q) (uj + iU-e(q)f 



6 2 + 



e -Pe{ q) + 



Here 



+ 47T 

92(6) = 4n 2 f(6) 



2 + 1T 2 )(9 2 + 4lT 
2 



2 1 



d9 e-^<V 

d9 e-? A < e+e ^ g 2 (9), 

3 

47T 2 g 47T 2 fl 2 " 

sinh 9 9 2 + n 2 



(158) 



(A0«(0 + 9 q ) + p^ c (9 + e q j) f(0) + 
1 c(9 + 9 q ) 



Ac(9 + 9 q ) df{9) 



d9 



(159) 



9 



I 2 + TI 2 )(6 2 + All 2 )' 

'vq s 



arcsinh 



A 



(160) 
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The next step in the low-temperature expansion is to subtract the contributions due to 
the partition function following (17), (11). The relevant contribution to the partition 
function is 



lim 35(h) 



k->0 



(161) 



Using that Ef ~ Ef + Ef 2 — Z\E^ X we obtain 

v 1 



+ 



e(g) + zO - e(q) 
v A 

(u; + zO-e(g)) 2 



d9 e-P* < e+e ^ g 2 (6) 



-PAe(q) 



+ 



+ OT n (w,g). 



(162) 



u; + z'O + e(g) 

Starting from Cfi{u),q) we can determine the contribution JF* in an analogous way. 
This then leads to the following result for the leading finite temperature contribution 
to the expansion (23) for frequencies u ~ A 

2v 



Cf(u,q) 



+ 



(oo + iO) 2 -e 2 (q) 
2vA([u + i0) 2 + e 2 (g)) /•« 



^ e -/3A c(e+e 9 ) ^ 



— zu 



— V 



e(q)((uj + z0) 2 -e 2 (q)) 

d9 e~^ K SM (a(cu, q, 9), 9 (u, q, 9), 9) 



. 2tt S(lu, q, 9)^/§ 2 (uj,q,9)-AA 2 
de e~ A ^ K SM (a(u;, q, 9), 9 (u, q, 9) 
Tl 2tt s(u, q, 9) V4A 2 -§ 2 {u,q,9) 



+ v 



de 



de. 



Ks M (e-,e + (q,e,e.),i 
[n(9,u) -u( q ,e,e_)}u{q,e,e. 



(163) 
t > 4A 2 

and 4A 2 cos 2 7 < s 2 (u, q, 9) < 4 A 2 respectively, and S is the contour from —00 — ry to 
00 — 27 parallel to the real axis. 



Here S + and Tj are the segments of the real axis characterized by s (cj,q 



5.3. Finite volume regularization 

A second way of regularizing infinities in matrix elements is to work in a large, finite 
volume R. As was pointed out in Ref. [25], up to corrections that are exponentially 
small in system size the functional form of matrix elements remains the same as in the 
thermodynamic limit. The main effect of the finite volume is to quantize to momenta 
or equivalently the rapidities, that parametrize the basis states. The quantization 
conditions in the finite volume are 

6*****% . . . 9 n ) bl .., n = S a b ^- 1 (9 nl ) . . . S^paOlflx . . . 9 n ) ai ... an , (164) 
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where 9ki = 9 k — 9\. In the one-particle sector we simply have 

J-R- sinl 

e « 

This is readily solved 



9j = arcsinh f j^^ ^j , j 62. (166) 

The density of such one particle states is 

A R 

Pl (6) = cosh#. (167) 

v 

In the two-particle sector the finite volume quantization conditions are obtained 
from the S-matrix eigenvalues S a (9) corresponding to spin singlet, triplet and quintet 
representations 

e iR± sinhfl! = e -iR± sinh0 2 = g(a) {Q 2 - d X ) , a = 0, 1, 2. (168) 

We have 

S^) = «±™ , (169) 

9 + Til 

In practice it is useful to consider the logarithmic form of the quantization conditions 

DA 

Y} S \9 U 9 2 ) = — sinh^- -J2 6iS) (*i - **) = 27rJ f } » ( 172 ) 

U Mi 

where S = 0, 1, 2, jj 5) GZ+^and 

2^ 



2 arctan 

5 (1) (6») = 2 arctan (J?-^ - 2 arctan (^j 



5^(9) = -2 arctan ( - ) . (173) 



7T / 

For later use we define the density of two-particle Bethe ansatz states with total spin S 

p?He h e 2 ) = iet dY 'lf' e ' ) . (174) 

o9 k 

In Bethe ansatz solvable models each solution of the quantization conditions (Bethe 
ansatz equations) gives rise to the highest weight state of an entire multiplet of the global 
symmetry algebra [58]. For the sigma model this means that each solution of the of (172) 
gives the highest-weight state of a 0(3) multiplet. The entire multiplet is constructed 
from the highest weight state by acting with the spin lowering operator. This leaves 
the spatial part of the wave function unchanged, and all states in the multiplet are 
characterized by the same set of quantized rapidities. 
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Recalling the action of spin operators on single-particle states (131), we may 
construct the following basis of two-particle eigenstates with definite values of total 
spin 

e 1 , e 2 - 2, ±2) = ^{iWn - 1 W 22 ± i|0i0 2 )i2 ± i|0i0 2 ) 2 i}, 

9 X , 9 2 , 2, ±1) = |{|^i^ 2 >i3 + IWai ± i\9i6 2 ) 32 ± i\e x e 2 )^, 
Oi, 2 ; 2, 0) = -L{2|^ 2 > 33 - IWii - IW22}, (175) 



9 h 9 2 ; 1, ±1) = ^{1^1^2)13 - I W31 ± ^1^2)23 T i^ra} 
9i, 2 ; 1, 0) = —j= 1 1 9\ 9 2 ) 12 — |0i02)2i}, 



(176) 



\9 X , 2 ; 0, 0) = -L{\9 X 9 2 ) 33 + \9 t 9 2 ) u + |6» 1 6> 2 ) 22 \. 



(177) 



The form factors in the basis (175), (176), (177) are readily obtained. Using the crossing 
relations for different rapidities 0i >2 ^ 9 3 

blb2 (9 x , 2 |$ a (O, 0)|6> 3 >6 = (0|$ a (0, 0)|6»! + in, 9 2 + in, 9 3 ) blb2b , (178) 

we obtain the following result for the form factor squares involving the two-particle 
singlet, triplet and quintet states 

r6 



J2 I Pi, 02; 0, 0|$ a (0, 0)\9 3 ) b \ 2 = - ip(9 12 )ip(-9 12 )(9 2 12 + n 2 ) 
b 6 

x ip(9 13 + m) 2 ^{9 23 + in) 2 , 



(179) 



E E 1(01,02; i^|$ a (o,o)|0 3 ) t 



b cr=-l 



7T 



- T ^(01 2 )^(02l)(013 + 02 3 ) 2 
X ^(013 + ^) 2 ^(0 23 + ZTr) 2 , 



;iso) 



E E l(0i,02;2,a|<|. a (o,o)|03) i 



6 <r=-2 



57T f 

12 



^(012)^(02l)(01 2 2 + 47T 2 



x^(^3+^vr) 2 7/;(0 2 3 + 27r) 2 . (181) 



In order to proceed we now assume that the form factors in a large finite volume have the 
same functional form as in the infinite volume up to exponentially small corrections. This 
is true for the quantum Ising model, for which the exact finite volume form factors are 
known, and support in favour of this hypothesis for general massive integrable QFTs has 
been provided in Ref. [25]. We are now in a position to evaluate the leading contributions 
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Co, C\ (22) to the low-temperature expansion (23) of the two-point function (134) of 
the nonlinear sigma model in the finite volume regularization scheme. The leading 
contribution is 



C?(u t q)KE*(u,q) + F*(u,q) 



1 1 



e(q) [id + iO — e(q) lu + iO + s(q) 

(182) 



(tu + tO) 2 -e 2 ( q y 

where the momentum is quantized q = 2nj/R with j an integer. The first subleading 
contribution is 

C?(u, q) = £? (u, q) + q). (183) 

Here S^(u, q) « E± + E± 2 (uj, q) — Z R E 0y i(u, q), where 

2 r R 



E« 2 (u,q) = J2E? s 2 (u,q) = J2[dx e**E{ 



5=0 5=0 
5 



2 b^s7?e^e 2 PM A 10 (01,92) 

x |(0 1 ,0 2 ;^ ( 7|$ a (O,O)|e3) 6 | 2 , (184) 

-!3Ac{e 3 ) -ixA[s{e 3 )-s{e 1 )- s {e 2 )} 

W(u, x, 0!, 9 2 , 9 3 ) = — Af ^TT , (185) 

u; + «d - A[c(6'i) + c(6> 2 ) - c(# 3 )J 

^ = 3^exp(-/?,(^)). (186) 

Here the sums are over solutions to the Bethe ansatz equations (166), (172). The 
contribution ^^(lj, q) is obtained from £^(u, q) using the symmetry 

F?{u>,q)=[£?(-L,,-q)Y. (187) 



5.4- Evaluating E^ 2 as R — > oo: Comparing the Finite and Infinite Volume 
Regularization Scheme 

In this section we compute E± 2 (u, q) as R — > oo. We will show that lim^oo E^ 2 (lj, q) — 
Eq ^uj, q)Zi is the same as £1,2(0;, g) — -Eo,i(o>, 9)^1 as computed in the infinite volume 
scheme. This then establishes again (as with the quantum Ising model) that the two 
regularization schemes are equivalent. However in this case we will have shown this to be 
true in a non-trivial setting: that of an interacting theory with non-diagonal scattering. 

To evaluate E^ 2 (uj,q) we first examine the contributions of each spin sector, i.e. 
the terms in the sum E± 2 (u,q) = J2s ^12 i^i <?) ( see Eqn. (184)). We will take the 
Fourier transform in q in the end and consider Ei 2 (u, x) for the time being. The latter 
takes the form 

E^,x) = IJ2Y. ^f" ! (|)o»th ! (|)cotf(|) , (188) 

z to to^ 2 pilhm (»u»2) 111 
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where 



g^{u, x, {6j}) = - tff(e 12 )tff{e 21 )tf>{e w + mf^e^ + ^) 2 

x coth 2 (^)tanh 2 (^)tanh 2 (^)^(^,g,{^}), (189) 



uWfax^Oj}) = ^-(9 2 12 + ir 2 )W(u,x,{0 3 }) 

7T 6 



f/«(o;,x,{^}) = -{9 lz ^02zfW{u,x,{9 j }) , 
U ( - 2 \u,x,{9 j }) = — (9 2 2 + Att 2 )W(u;,x,{0 j }). 



(190) 



Like with the Ising model in Section 3.6, our aim is in the R — > oo limit to convert these 
sums into integrals. To do so we must first isolate and subtract out the singular terms 
in E^^ou, x) in order to define finite integrals. Expanding Ef-^iu, x) in 9 3 about 9\ and 
02, we obtain its singular pieces 

= ^E E , n , mr n 7T E r l s) (^^W) 



2 V^Plfa) pf^lA)^ 

y E {a) 



(191) 



a=l 



rfV, 



16 t/( g >(^,s,fli,fl 2 ,fl 2 ) cosh(fl 2 ) cosh(fl 3 ) 
^ (0 2 i2 + * 2 )(Q 2 i2 + 4vr 2 ) (sinh(# 2 ) - sinh(# 3 )) 2 

+ (0i ~ 2 ) , 

16 1 cosh(fl 3 ) 



(192) 



7T 



4 % 



+ 7T 



-2UW(u,x,6 1 ,6 2 ,6 !i ) 

+ (^2 + ^21 J- 



j -r 4vr 2 )6>2i sinh(6> 3 ) - sinh(6' 2 ) 

#21 



7T 



sinh(# 2 ij 



7T 



6»3=e2 



+ 



(193) 



The singular pieces, as indicated by the introduction of the quantities E^ g {u,x) come 
in the form of both single (a = 1) and double poles (a = 2) in # 3 about 9\ and 9 2 . To 
regularize the sum we add and subtract the singular pieces from E^ 2 {uj,x\. 

'GW(w, x, {6j}) tanh 2 (^) coth 2 (^) coth 2 (^f) 



_ 771 H£> 



e 3 e^e 2 



HS 



Pl (e 3 )p { 2 s \e 1 ,e 

t[ S) (uj,x,9 u 9 2i 9 3 



-E 



— 1 Pl (e 3 ) P2 s \e 1 ,e 2) 



singl, 



Ju, x) 



2 \lo,x) + E^ n s lgl 2 (u, x). 



(194) 

The first term in the above equation is singularity free (i.e. the summand is finite as #3 
approaches either 9\ or 9 2 ). We can then take R — > 00, turning the sum into a principal 



Finite Temperature Dynamical Correlations in Massive Integrable QFTs 40 

value integral after which the integration contours can be modified so that they deform 
about the singularities in the same fashion as was done for the Ising model (see Section 
3.6). The result of doing so is 

*ffiu>, x) = / dJ ^G'^, *, m trf(^) coth*(^) caU.'(*!±i2) 

+ / '^^1 i ( e 32 )i(93i) G< S W,{%» • (195) 

We note that the principal part integrals over the single-pole terms vanish. If we now 
sum over the spin sectors, S, then take the limit R — > oo and finally Fourier transform 
with respect to x we find 

+ 87T 2 / M )f*®* w(u>, q, 9 U 9 2 , 9 3 )5(9 12 )6(9 13 ), (196) 



(2tt) 3 

where #2, #3) is defined in (153) and W(u, q, 9i, 9 2 , 9 3 ) is 

vCr< n n /n e-^^2n5(vg + A(s(9 3 ) - 3(9,) - s(9 2 ))) 

W{u, q, 0i, 9 2 , 9 3 ) = — — . 197 

+ id - A[c(9i) + c(9 2 ) - c(9 3 )\ 

We see that we have reproduced El° 2 n (uj, q) as given in Eqn. (152) plus a disconnected 
term. The remaining disconnected terms together with a term proportional to the 
partition function are found E^ gl 2 (a>, q) which we now turn to evaluate. We are able 
to carry out the sum over 9 3 courtesy of the identities 

i + s (s) (M 



y-v COsll(# 



E 



Pl (e 3 )(sinh(e 3 )-siiiii(e 2 )) 

cosh(6 ) 2) cosh(6> 3 ) 



Pl (03)(sinh(0 3 )-sinh(02)) 2 

03 



AR cosh(# 2 



- (198) 
2v 1-Re S( s )(9 21 ) ' 1 J 



Here 9\ ^ 9 2 £ S indicates that that 9i )2 are solutions of the Bethe ansatz equations (172) 
in the spin-S sector and S^ s \9) are the S-matrices (169). The first identity is established 
by replacing 9 3 by its corresponding integer through the quantization conditions (166) 
and carrying out the resulting sum using 

1 2z °° 1 

»'W = -+ 7 E^ («») 

k=l 

Finally, the Bethe ansatz equations (172) in the two-particle sector are used to rewrite 
the result. The second identity in (198) can be established by using the derivative of 
the first. 

The second identity in (198) allows us to evaluate the double pole term in 
E^ gl 2 (u, x) with the result 

2S+1 ^ ARcosh(9 2 )W(uj,x,9 1 ,9 2 ,9 2 ) 



4 2 >,*) = ^ E 



3 " »,fe S pf\0U02) 
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_2S+1 &Rcosh(9 2 )W(uj,x,9 1 ,9 2 ,9 2 ) 



3v e^es pi\ex,e 2 ) 
2S + 1 ^ w{uj,x,e 1 ,e 1 ,e 1 ) 



_ 25+1^ 



3 ^ Pi(fli) 

"^~3-? M^O • (2 ° 0) 

Here we need to subtract the term 9\ = 9 2 only in the triplet sector because in the 
S = 0,2 sectors solutions of the Bethe ansatz equations (168), (169) with coinciding 
rapidities do not occur. In the triplet sector the spatial part of the wave function has 
to be antisymmetric (as the spin part is and we are dealing with bosons), which forbids 
solutions with coinciding rapidities. Turning the above integral back into a sum over 
rapidities subject to free quantization conditions (166) we arrive at 



p(2) (,, T \ = 2S+1 ST exp (ix^s(9 1 )) ^ 



-/3Ac(0 2 ) 



3 ^p 1 (9 1 )[u + i6-Ac(9 1 )} 
-hi [ ^ L W(u,x,9 1 ,9 1 ,9 1 ). (201) 

J 27T 

Carrying out the Fourier transform in x and summing over spin sectors we obtain 

£ & ( £ g (^, q) = O, q)Z* - J ^ W(u, q, 9,, 9 U 9,), (202) 

where W is the spatial Fourier transform of W. We thus obtain a term proportional 
to R as R — > oo but which will be cancelled off by a corresponding term, Eq^u, q)Z^, 
arising from the expansion of the partition function. 

The first identity in (198) allows us to evaluate the single pole term in E^ gl2 (u, q). 
Unlike the double pole term, this leads to an expression entirely finite in the R — ► oo 
limit: 

2 /•dMM^s W(u,x,9 h 9 2 ,9 3 ) 



5> s S g (^,z) = 167r 2 ! 

c J 
l(ff 2 + 

sinh(0 2 ij 



(2VT)3 (^ 2 + 7r 2 )2( ^ 2+47r2) 

'^(032) {2^(2 ^ -fp + - tt 2 )) + (0 2 2 + 7 r 2 )(2 7 r 2 + 2 2 )} 



+ 7T 2 21 (0 2 2 + 7r 2 )5 / (032 

This then allows us to write down a complete expression for E^ 2 (uj, q) as R — > 00: 
E&(u> t q) = lJ d -^^ W(u,q,9 1 ,9 2 ,9 3 )L(9 1 ,9 2 ,9 3 ) 



(203) 
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We see that this expression agrees with that of Ei 2 (uj,q) using the infinite volume 
scheme, see (152), (153) and (155). 

5.5. Resummation 

In the above we have calculated the first two terms in the expansion of the dynamical 
susceptibility 

X, (u, q) = C V, q) + C?(u, q) + .... (205) 

The calculation of Cf is not possible in practice because it involves the 5-particle form 
factor, which is tremendously complicated for the nonlinear sigma model [7]. We have 
shown that Cf(u,q) exhibits a quadratic divergence when the frequency approaches 
the mass shell lo — > e(q), whereas C*(dJ,g) is only linearly divergent. We have argued 
that the higher order (in exp(— A/T)) terms in (205) exhibit stronger divergences and 
hence a resummation is required to get meaningful results for uo m e(q). Following the 
resummation procedure set out in section 2.1 we determine the quantity 

£ (I W = -g^, (206) 

and then use it to obtain a resummed low-temperature approximation (see (27)) 

,(Dr„ a\ = ff^g) 



X^,q) = - ^°^ 1U ., ~ , (207) 



SM(u,q) = - I-J-^I m x«(u,,g). (208) 

7T 1 — e T 

6. Results for the Low- Temperature Dynamical Susceptibility of the 0(3) 
Nonlinear Sigma Model 

The leading order result S^'(uj,q) is most easily calculated using the infinite- volume 
regularization scheme and carrying out the integrals in (163) numerically. Evaluation 
of S^(u), q) for low temperatures shows that, as expected, the T = delta function at 
uj = e(q) broadens with temperature. We find that the resulting peak scales as 

A /A\ 
peak height oc — ex P (-J, 

peak width oc — exp(-— V (209) 



A V T 

In order to exhibit the evolution of the structure factor as a function of frequency for 
fixed momentum q with temperature it is therefore useful to rescale both the frequency 
axis and the structure factor. The result is shown in Fig.9 for a temperature range 
0.1 A < T < 0.3A, which corresponds to approximately a factor of 2000 difference in 
peak height (and width). We see that the lineshape is asymmetric in frequency, with 
more spectral weight appearing at higher frequencies. The asymmetry increases with 
temperature. This effect is most easily quantified by comparison with a Lorentzian 
lineshape, which is done below. 
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Figure 9. Rescaled structure factor for the 0(3) nonlinear sigma model for three 
different temperatures. 



At very low temperatures we find that the line shape at q = is well-approximated 
by a Lorentzian 

^W^ C-WtPgT (210) 

A comparison of S^(u,0) to (210) is shown in Fig. 10 for two temperatures. We see 
that the agreement improves with decreasing temperature. 

6.1. Comparison to Semiclassical Results 

In Ref. [39] Damle and Sachdev carried out a semiclassical analysis of the dynamical 
structure factor of the 0(3) nonlinear sigma model. In contrast to the case of 
the quantum Ising model, the scattering matrix is not taken into account fully but 
approximated by its zero rapidity limit 

S${8) — SadSbc. (211) 

We note that compared to the Ising case this additional approximation may impose 
tighter restrictions on the window of applicability of the semiclassical result. Damle 
and Sachdev find the following form for the dynamical structure factor 

r 00 tit 

S sc (u, q) = / — dx e^-^Kix, t)R(x, t) , (212) 
J~oo 2tt 

where the relaxation function R(x, t) is determined numerically and 

K(x,t) = -K (\A\^/(x/vy-t*) . (213) 

7T 



Finite Temperature Dynamical Correlations in Massive Integrable QFTs 



44 



1 1 1 1 1 1 1 1 1 


1,1,1, 


— T=0.08 A 

— T=0.1A 

i 





1.005 



0.995 



_l i I i I i I i L_ 

-10 -5 5 10 

(co-A)/r 



Figure 10. Comparison of the structure factor to the Lorentzian approximation (210). 



Here Z is a normalization factor. At sufficiently low temperatures and when \vq\ < VTA 
the semiclassical result is approximately Lorentzian in form [39, 38] 

vZ ct/r Q 



S sc (u,q) « S Lov (u,q) 



ne(q) (uj-e(q)) 2 + (a/r y 



where a ~ 0.72 [38] and 



To 



7T 



A/T 



3T 



(214) 



(215) 



We have shown above that our result is well approximated by a Lorentzian at low 
temperatures. However, unlike (214), the best fit to our result involves a temperature 
dependent gap, see (210). The width of the Lorentzian is quite close to Damle and 
Sachdev's result, e.g. we find that T(O.IT) = 0.736/r^,. However, our value of a is 
found to increase as T — > 0. In Ref. [41] an analytic expression for the relaxation 
function R(x, t) was given 

d<fi (1 — 7 2 ) cos(x sin 0) 



R(x,t) = C 



2n 7 2 + 27 cos 0+1 



x exp 



COS0)|t 



-, /e 



where C is a normalization constant, 7 



3' x 



u erf(u] 



x lioi t = t/r , u = x/t and 



£0-3 



2nv 2 
~k^T 



3 A/T 



(216) 



(217) 
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At low temperatures and for u sufficiently close to e(q) the Fourier transform (212) is 
well approximated by [39] 

dt 



S sc (u, q) 



2tt 



dx e iuJt - iqx K(x,t)R(0,t) 



Carrying out the integrals we obtain 



S sc (uJ,q) 



ZCv 



where 



r(T) 



7re(q) 
7+; + 2 



Re 



T(T)+i(u-e{g)) 
(lj - e{q)Y + r 2 (T) 



1 + 



7 



27^ 



a = 1 - i^r (oJ - e(q)) . 



(218) 



(219) 



(220) 



While this expression is very close to being a Lorentzian sufficiently far away from the 
mass shell, it exhibits a square root divergence as uj — > e(q). This precludes a comparison 
with the results of our low-temperature expansion. 



7. Summary and Discussion 

In this work we have proposed a general method for determining frequency and 
momentum dependent two-point functions of local operators in massive integrable 
quantum field theories at low temperatures. We have applied this method to 
the calculation of response functions in the disordered phase of the quantum Ising 
model and the 0(3) nonlinear sigma model. The methodology in this particular 
application possesses two crucial ingredients. First we have shown that there exists a 
systematic expansion of the spin-spin response function in terms of the small parameter 
exp(— A/T), where A is the spectral gap and T is the temperature. While such 
expansions have been known to exist previously [23], they were limited to certain 
correlation functions. Here we have shown that the spin-spin response function can be 
described sensibly by performing a low temperature expansion of the spin-spin response's 
"self-energy" (or a quantity analogous to such). The second ingredient was a procedure 
to make sense of infinities that appear in the squares of matrix elements that arise for a 
Lehmann expansion of the correlators. In this paper we have developed a new regulator 
for the infinities that appear when working in an infinite volume and shown that this 
regulator reproduces results found by working in a large, finite volume and taking the 
thermodynamic limit only at the end of the calculation. We have accomplished this 
both in the quantum Ising model, a free fermionic theory, and the 0(3) NLSM where 
the elementary excitations are strongly interacting with a non-diagonal, momentum 
dependent scattering matrix. 

A number of open problems remain on the technical level. In the present manuscript 
we have not analyzed the situation where we need more than one k parameter in the 
infinite volume regularization scheme (see e.g. (59)). It is important to establish the 
equivalence of infinite and finite volume regularization schemes in this more general case 
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as well. We also have not presented a general proof that all terms (17), (18) in the low- 
temperature expansion (19) are finite. We hope that these questions will be addressed 
in future work. 

At zero temperature the response functions in both the Ising model and the 0(3) 
NLSM are dominated by a delta function peak at the position of the single particle 
dispersion. We have determined how this peak broadens at finite temperatures smaller 
than the gap. Our main result is that the lineshape at T > exhibits a pronounced 
asymmetry in energy that increases with temperature. At very low temperatures 
T < 0.1 A our results essentially reduce to those of previous semiclassical analyses 
[28, 39, 38], which concluded that the lineshape is to a good approximation Lorentzian. 
This shows that while the semiclassical approximation gives a good account of the width 
and height of the lineshape at T > ||, it is quantitatively accurate with regards to the 
overall lineshape only at very low temperatures. 

A second important feature seen in our low-temperature expansion is what is known 
as the temperature dependent gap. For sufficiently high temperatures the maximum of 
the lineshape is seen to shift upwards in energy when compared to the T = gap. 
We found that this phenomenon emerges in the quantum Ising model once subleading 
terms in our low-temperature expansion are taken into account. The calculation of these 
terms involves five-particle form factors. While these are known for the 0(3) NLSM as 
well, their complexity puts the calculation of subleading terms in the low temperature 
expansion beyond the scope of this work. 

The 0(3) NLSM describes the scaling limit of integer spin Heisenberg chains. 
While the agreement between the field theory and the Heisenberg lattice model is 
best for large spins (and low energies), the sigma model has been found to provide 
a reasonable approximation of the structure factor even in the extreme 5 = 1 case [59]. 
A number of inelastic neutron scattering experiments have measured the temperature 
dependence of the dynamical structure factor of quasi one dimensional spin-1 Heisenberg 
magnets[34, 35]. Our finding of an asymmetric lime shape are relevant to these 
experiments. In particular, the excess of spectral weight at high energies reported in [34] 
for CsNiCl3 should at least be partly accounted for by a lineshape asymmetric in energy. 
However, a quantitative comparison of our theory to inelastic neutron scattering data 
on CsNiCls is precluded by the presence of a Neel transition, driven by non-negligible 
interchain coupling, at a temperature comparable to the gap. It is likely that the 
interchain coupling will affect the precise lineshape in a significant way [60]. Similarly, 
in YBaNiOs a quantitative comparison to the finite temperature structure factor is not 
straightforward because of the presence of an exchange anisotropy. 

The presence of an asymmetric lineshape even at low temperatures is expected to be 
a general feature in quantum magnets that support coherent single particle excitations 
at zero temperature. Theoretical studies of the alternating Heisenberg chain [61] suggest 

|| Indeed, as was shown in Ref. [36], the width of the lineshape obtained in the semiclassical 
approximation is in agreement with experiments on a number of different materials (with gaps that are 
small compared to the magnetic exchange constant) even at elevated temperatures TkA. 
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that asymmetric lineshapes occur in dimer systems. This has been confirmed by recent 
experiments on the quasi one dimensional alternating Heisenberg chain copper nitrate 
[62]. It would be interesting to investigate to what extent the same holds for the two 
and three dimensional cases [63]. 

The methods we have developed have a wider scope for applications. With regards 
to low temperature dynamics in gapped integrable models, it would be interesting to 
investigate the case of the spin-1/2 Heisenberg-Ising chain. Here the dynamical structure 
factor has been measured for a number of materials [64]. In the limit of large gaps the 
transverse component of the structure factor has been determined by diagrammatic 
methods [65]. The excitation spectrum is quite different compared to the 0(3) NLSM 
and the disordered phase of the quantum Ising model in that the lowest excitations 
are two-parametric. Concomitantly the dynamical structure factor is dominated by an 
incoherent two-particle continuum at zero temperature. At finite temperature a low- 
frequency resonance, known as "Villain mode" develops [66]. It should be possible to 
adapt our methods to analyze this case. 

Our method of regularizing general matrix elements could be useful for studying 
non-integrable perturbations of integrable models, which is an area of considerable 
importance [67]. In fact shortly after our work a preprint by Takacs appeared 
(arXiv:0907.2109), which addresses this problem for the double sine-Gordon model and 
obtains an independent derivation of the finite volume regularization scheme. 
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Appendix A. Crossing Relations 

In this Appendix we summarize identities first given by Smirnov in [1] that that allow 
us to analytically continue form factors. Let A = . . . , 9 n } with 9\ < 9 2 < ■ ■ ■ < 9 n 
and B = . . . , f3 m } with /3i < (3 2 < ■ ■ ■ < (3 m and introduce notations 

Z[A] ai ...a n = Z ai (9i)Z a2 (9 2 ) ... Z an (9 n ) , 

Z^} an ... ai = Zl n {9 n )Zl_ x {9 n ^) . . . Zl{9,) . (A.l) 

Now let A\ and A 2 be a partition of A. As a consequence of the periodicity axiom we 
have 

z[l] ai ... an = z[T 2 ] Cl ... Cr z[Al] Cr+1 ... Cn s(1\aI):\:: c z, (a.2) 

where S( A \ Ai) is the product of two-particle scattering matrices needed to rearrange the 
order of Faddeev-Zamolodchikov operators in Z[A] to arrive at Z[A 2 ]Z[A-\\. Similarly 
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we have 

Zt[£] 6ro ... 6l = tf[%] dm ... da+1 tf\<B 2 ] ds ... dl S(%\Wct- ( A -3) 
Finally we define 

n 

8[A,^W- r = 6„.,„\\-l~6{0 J - 6 (A.4) 

o 1 ■ ■ ■ "m A A 

For a local operator (9 we then can analytically continue form factors as 
(0\Z[1] ai ... an 0(0,0) Zt[% m ... bl |0> 



C1 ...c r 
d s ...d„ 



A=A 1 UA 2 
B=B 1 UB 2 



x (0|Z[X + zO] Cr+ ,.. Cn O(0,0)Zt[ < B 1 ] dm ... d3+1 |0). (A.5) 
Similarly, we could choose to analytically continue to the lower half-plane 
(0\Z[1] ai ... an 0(0,0) Zt[% m ... 6l |0> 

A=AjUA 2 
B=B 1 UB 2 

x (OlZ^x -iO] Cr+1 ... Cn O(0,0)Zt[^ 1 ] dm ... d . +1 |0). (A.6) 
For semi-local operators the above identities needs to be modified as discussed in [1]. 

Appendix B. Products of Form Factors for the Ising Model 

Let us consider the product of form factors 

K0 3 KO,O)|0 2 0i>| 2 = ]hn(9 3 \a(0,0)\e 2 e 1 )(e 1 ,9 2 \(T(0,0)\9 3 + K). (B.l) 

re— >0 

Using the crossing relations (57) this can be rewritten as 

| (0 3 |a(O, 0) |0 2 0!) | 2 = lim{ ((0 3 |<7(O, 0) \9 2 - iO, 9 X + iO) - 27ra[5(0 32 ) + 5(0 31 )]) 

X ((0! + zO, 02 - lO|a(0, O)|0 3 + K) + 27Tff [5(032 + «) + 5(031 + «)]).} 

(B.2) 

Multiplying out the various terms then gives 

1(03^(0, 0)|0 2 0i)| 2 = lim{ (0! + z'0, 2 - zO|<t(0, O)|0 3 + k)(9 1 - z'0, 2 + iO\a{0, 0)|0 3 )* 

+ 2tto- [5(032 + k) + 5(03i + k)] (0i - zO, 2 + z'0|a(0, 0) |0 3 >* 

- 2na [5(0 32 ) + 5(0 31 )] (0! + zO, 2 - z0|tr(0, 0) |0 3 + «> 

- (27T.T) 2 [5 (032 + «) + 5(031 + «)] [5(0 32 ) + 5(0 3 l)]} 

= lim [r conn + r^ 8 ' 1 + r dis ' 2 ] . (b.3) 

The connected part 

pconn = (0 1 + iO ,0 2 -iO|(7(O,O)|03 + «) 

x (0!-zO, 2 + z'O|<7(O,O)|0 3 )* (B.4) 
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does not contain any divergent pieces and the limit k — > can be taken straightforwardly. 
The product of delta-functions gives 

(B.5) 
(B.6) 



+ S(9 31 )S{9 32 + k)+S{9 32 )S(9 31 + k) 
where 9jk = 9j — 9k- The cross-terms are 



^dis,2 



2,710' r 



cross, 1 



27iaT 



cross, 2 



where 



r c ross,i = 5(032 + k)(9i - if]i, 9 2 + ir) 2 \a\9 3 )* 

- 5(032) (0i + iVi, 02 - ^ 2 k|0a + re), 

Tcross.2 = 5(031 + «) (01 - «7l, 02 + Mfcl^)* 

- 5(0 3 i)(0i + iv'i, 02 - ^ 2 k|03 + re). 

Here r\\ i2 are positive infinitesimals and we are interested in the limit 771,2 — » at fixed 
k. Using the explicit form of the form factors we obtain 



(B.7) 



r cr oss,i = -icr5(9 32 + k) tanhf 



0i2 + irji + ir}2 



coth 



0i2 + re + irji \ 
9. J 



coth 



K - vr} 2 



+ io~5(9 32 ) tanh 



012 + iv'i + iv'2 



coth 



012 ^ 4" irji 



coth 



K + irj 2 \ 



)■ (B-8) 



Using that 



tanh 
~ 1 



0i2 + irji + ir l2 



sinh 9 



12 



2 

(re - ^ 2 
this can be simplified to 

r cr oss,i = - 2iaS(9 32 + k 

+ 2ia5(9 32 ) 



coth 



012 + re + z?7i \ 

2 / 
1 1 



012 012 + re + 2771 
1 1 



(B.9) 



1 



1 



k - ir) 2 sinh 9\ 2 9 12 9 12 + k + 2771 . 

1 1 1 

+ — — — + 



. k -\- irj 2 sinh 9\ 2 9± 2 9 12 - k + irji J 
2na5(K) [5(9 32 + k) + 5(0 32 )] + ^3^3 P 



2za [5(0 32 + re) + 5(9 32 )} 

S(9 32 ) 5(9 32 + k 



K z + 775 
1 1 

sinh 9 12 9 12 



+ 2ia 

In the limit k — > this becomes 

1 cross, 1 -> 47ra5(^32)5(«) + 2ia 



012 - re + 2771 i2 + k + ir]i 
S(9 32 ) 



(B.10) 



- 2iaS'(9 32 ) + m6(9 32 ] 



S(9 32 + ac) 
6»i2 - k + id 6»i2 + « + iO 
1 1 



sinh 0i2 



(B.11) 
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6(0 



31; 



5(0 3 1 + «) 



2^'(03l)+4i^(e3i; 



012 + « + ^0 i2 - « + i0 

1 1 



(B.12) 



_sinh^i2 0i2. 

We note the emergence of a term involving the derivative of the delta function. Our 
final result for r dls ' 2 is then 

lim r diSi2 = limjsvrV [5(9 31 ) + 5(0 32 )] 5(k) 



+ Aitia 2 [5'(0 31 ) ~ S'(0 32 ] 
+ 8ma 2 [5(0 3 i) + 5(0 32 )] 



1 



+ 4nia 



sinh 9 12 0i 2 

S(9 3 l) + 5(032 + «) , 5(031 + «) + 5^ 



6»i2 + « + zO 6>i2 - « + i0 

The infinite volume regularization for the form factor squared is then 



}. (B.13) 



1(03^(0, O)|0 2 0i)| 2 = cr 2 tanh 2 f^) coth 2 (^i±^ coth 2 ^ 23 



iO 



2 J V 2 
+ lim [r^ 8 ' 1 + T dis ' 2 ] , (B.14) 

where r dlSi1 and r dls ' 2 are given by (B.5) and (B.13) respectively. For our purposes we 
only need the part of r dls = r^' 1 + r dls ' 2 symmetric under 9\ <-> 9 2 as the remainder of 
the integrands we consider all have this symmetry. The symmetric part is 

r dlS (6»i,6» 2 ,0 3 ) + r dlS (6> 2 ,01,0 3 ) „ 2-2r/ n [f/fl \.X/fl M 

- An a d{K) [d(9 31 ) + d(9 32 )\ 

+ 87rV 5(0 3 i) 5(9 32 ). (B.15) 



lim 



Appendix C. Products of Form Factors for the 0(3) nonlinear a-Model 

In this Appendix we determine the form factor squared 

r(0i,0 2 ,0 3 )= U(0i,02|$ a (o)|0 3 > 6 | 2 . (c.i) 

Using (150), which is a consequence of the crossing relations discussed in Appendix A, 
we obtain the decomposition 

r(#i, 2 , 9 3 ^ — Ktti ^ conn (()_ /3_ Q„\ _L r dis ^ id, /3_ OA _L r dis > 2 C 



Here 



= lim [r™(0i, 2 , 3 ) + r* 8 - 1 ^, 2 , 3 ) + r dis > 2 (#i, e 2 , 0, 



r™= fel 6 2 (0i-^i,02 + ^2|$ a (o)|0 3 ) 6 * 

X M 2 (01+^l,02-^2|$ a (O)|0 3 + K) 6 , 
T^' 1 = ^(2VT) 2 S a >2l) [5(032)5(032 + «) + 5(031)5(031 + «)] 

+ ^(2VT) 2 S#(0 21 ) [5(031)5(032 + «) + 5(0 32 )5(031 + re)] , 



(C.2) 



(C.3) 



(C.4) 
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b 

+ ^27r5(0 3 i + «) ba(Oi -iVi,02 +im\* a (O)\0i ~ «)2 

6 

+ ^27n5(0 32 ) b a(# 2 - i7te,0i + zr/i |$ a (0) |0 2 + n) b 

b 

+ ^2^(^31) ab(0 2 -iV2,9i +iVi\* a (O)\0i + ■ 

b 

Using crossing the connected contribution can be written as 

peon* = J- (0|$ a (0)|^ + for" 2 + *7T + , fo)^ 
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(C.5) 



to 



X (0|$ a (0)|^ + 27T+, 2 + ITT", 3 + 



(C.6) 



Using the explicit expression (138) for the three-particle form factor this can be brought 
in the form 



lim r conn 

K,T]1,2— »0 



7T 



0? 2 + ,r 2 



tanh 



4/^12 



01 



(^13 + i0) ; 



(023 -^0) 2 ,,4/ 

COth (—2— J (0| + vr 2 ) 2 C ° th I 



2 

023 - *0 



(C.7) 



(0?3 + 

The connected r conn contribution is now in a form suitable for further analysis. For our 
purposes it is sufficient to determine the parts of p 118 ' 1 , r dls ' 2 symmetric in 6\ and 9 2 

r dis ' j (0i,02,0 3 ) + r dis 'j(0 2 ,0 1 ,0 3 ) 



^dis,j , 



h, 02, 03 



lim 



We therefore concentrate on these parts only in the following. We have 
limrf' 1 = lim V4tt 2 S(k) ReS$(9 21 ) [S(9 32 ) + S(9 31 )} 



j = 1,2. (C.8) 



24vr 2 5(9 31 )S(9 



32 , 



(C.9) 



Using the explicit forms of the three-particle form factors and proceeding along the same 
lines as for the Ising model, we obtain after some lengthy calculations 

2rf 2 (0 1 ,0 2 ,0 3 ) = ^(27t) 2 5(0 32 )5(k) Re [l - ^ 6 6 (0 21 ) 



+ ^(27r) 2 «5(0 31 )5( K ) Re [l - S, 

b 

+ ^4tt5'(0 32 ) lm^(0 12 ) 

b 

+ ]T(47r) 2 cK0 32 )(l-<U) Re 

b 

+ $>vr5(0 32 ) ImS a t(0 12 ) 



abfn 



1 



021 + 



'21 



2m 
20 



21 



Sinh 021 021 0?! + 7T 2 



+ J](27r) 2 5(0 3 2 + k)5{9 31 ) Re [l - S a a b b (9 21 )} 
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+ ^(27r) 2 5(# 3 i + «)5(032) Re [1 - S, 



abfn 



^31 + »r 

#31 +« 



'12 J 



tab I 



'31 



7 2- 



(CIO) 



This is simplified further to 

4tt 2 [6{9 32 ) + 6{9 31 )}6{k) 



limrf' 2 - 



-i67r 4 [m 2 )-mi)]- 

+ 327r 2 5(^ 32 )5(^3i) 

<5(0 32 ) +5(^1) 



3 - 5>eS: 6 6 (0 21 ) 

#12 



'12 



+ 7r 2 )(fl 2 2 + 47r 2 ) 



+ 16tt" 



^2 + 



4vr 2 # 



4vr 2 # 12 



sinh 9 12 9\ 2 + 7r 2 



2/32 
12 



(c.n; 



(^2 + ^(^+4^) 

Our final result for the symmetrized (in 9\ and 9 2) disconnected parts of the form factor 
squared in the infinite volume scheme is 

12^ 2 [S(9 32 ) + S(9 31 )]S(k) 
+ 87r 2 5(9 32 )5(9 31 ) 
- 16tt 4 [S'{9 32 ) - S'(9 31 ) 



lim 



pdis,l _|_ pdis,2 



9 



12 



(0 2 2 + 7r 2 )(0 2 2 + 47r 2 ) 



+ I67H 



S(9 32 ) + S(9 31 ) 
h + ^)(9f 2 + A^) 



4n 2 9 



12 



4vr 2 ^ 2 



12 sinh 12 ^2 + 7T 2 



. (C.12) 
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